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I.  Introduction 

Near-infrared  light  offers  unique  possibilities  for  detecting  physiological  information  from  deep 
within  biological  tissues  without  injections,  surgery,  dangerous  radiation,  or  great  expense. 
Although  it  has  been  known  for  some  time  that  non-ionizing  near-infrared  (NIR)  light  penetrates 
a  few  centimeters  into  biological  tissues,  recent  modeling  advances  have  promised  hope  for 
clinical  applications.  NIR  light  provides  the  concentrations  and  oxygenation  state  of 
hemoglobin,  as  well  as  the  scattering  properties  of  the  tissue  and  thus  provide  criteria  to 
characterize  breast  lesions  in  vivo.  Optically  monitoring  physiological  changes  in  breast  tissue 
provides  a  new  way  to  look  at  both  disease  progression  and  prevention.  Our  prototype  NIR 
instrument  has  observed  changes  in  tissue  physiology  that  are  characteristic  of  disease 
progression  and  treatment.  We  have  also  observed  changes  in  breast  tissue  physiology  resulting 
from  hormone  replacement  therapy,  menopausal  status,  menstrual  cycle  stage  and 
chemotherapeutic  intervention.  Because  our  technique  is  safe,  non-invasive  and  comparable  in 
cost  to  ultrasound,  patients  can  be  studied  as  often  as  necessary  without  risk. 


II.  Body 

Il.a  Items  in  the  Original  Specific  Aims 
The  original  specific  aims  of  this  project  are: 

(1)  Optimize  the  instrumentation  and  theory  currently  used  to  extract  the  absolute  optical 
absorption  and  reduced  scattering  parameters. 

(2)  Perform  clinical  measurements  of  breast  tissue  optical  properties  in  order  to  determine 
exact  absorption  and  scattering  factors  in  different  tissue  regions. 

(3)  Compare  optical  property  measurements  to  gross  pathology,  histopathology,  ultrasound, 
and  mammography  results  in  order  to  determine  absorption  and  scattering  factors 
characteristic  of  tumor  and  normal  tissue. 

Items,  1  &  2  have  been  given  sufficient  attention  to  warrant  several  publications  over  the  cycle 
of  this  grant.  The  third  item  is  being  written  in  another  publication  as  a  detailed  case  report 
describing  the  monitoring  of  chemotherapy  in  a  breast  tumor. 


Il.b  Progress  on  Statement  of  Work 

Below  I  list  progress  on  each  of  the  items  listed  in  the  3-year  statement  of  work  as  originally 
described  by  Dr.  Joshua  Fishkin  (page  10  of  original  proposal).  After  taking  over  his  original 
grant  in  1999,  new  opportunities  arose  that  allowed  me  to  place  a  greater  emphasis  on  clinical 
measurements  as  opposed  to  developing  more  sophisticated  computational  models.  Given  the 
scarcity  of  clinical  data  available,  I  elected  to  pursue  this  option  at  the  expense  of  the 
computational  options. 

l.a  Selection  and  addition  of  new  laser  diodes,  modification  of  instrument  software. 
completed.  We  have  constructed  2  more  instruments  that  feature  10  laser  diodes 
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(upgraded  from  4).  The  new  systems  also  feature  a  method  providing  continuous 
coverage  from  650  to  1000  nm  (Paper  #4).  The  instrument  software  has  been  completely 
re-written  to  include  these  new  features. 

1  .b  Clinical  measurements:  completed.  Initial  clinical  studies  were  conducted  at  BLI.  Only 
non-invasive  measurements  were  performed  (see  clinical  summary). 

l.c  Accuracy  studies  of  one  and  two  chromophores  systems:  completed.  The  FDPM 
technique  reproduced  the  absorption  and  scattering  properties  of  a  tissue-like  water-based 
phantom.  For  higher  absorption  phantoms,  the  FDPM  technique  was  less  accurate. 
However,  the  SFDPM  technique  (Paper  #4)  greatly  enhanced  our  accuracy,  and  reduced 
errors  from  instrument  calibration. 

1. d  Spherical  &  cylindrical  phantom  studies',  not  completed.  Constructing  solid  phantoms 

proved  difficult.  A  company  (Medlite)  was  contracted  to  construct  these  phantoms. 

2.  a  Clinical  measurements:  completed.  See  clinical  summary  for  details. 

2.b  Studies  with  inhomogeneous  phantoms:  completed.  In  order  to  assess  the  accuracy  of  the 
physiological  measurements,  I  elected  to  compare  my  measurements  against  known 
normal  breast  physiology  (as  opposed  to  tissue  phantoms).  This  comparison  is  more 
useful  than  phantom  studies  because  a  phantom  can  never  reproduce  exact  human 
physiological  conditions  (see  clinical  summary). 

2.C  Accuracy  study  of  S  &  4  chromophore  systems:  completed.  Using  the  SSFDPM 
technique,  we  reproduced  the  spectra  of  a  4  chromophore  system.  Details  of  this  work 
were  reported  (Paper  #4). 

3  a.  Construction  of  tissue  simulating  phantoms:  not  completed.  Heterogeneous  phantoms 

were  not  constructed  (see  1  .d). 

3b.  Clinical  measurements:  completed.  A  clinical  system  now  resides  permanently  at 
UCIMC.  Our  clinical  systems  have  been  used  in  normal,  cancer,  and  chemotherapeutic 
monitoring  studies.  Results  of  these  measurements  have  been  the  subject  of  three  peer- 
reviewed  papers  (paper  #1-3),  with  more  planned.  The  results  of  these  measurements 
have  shown  that  the  FDPM  technique  can  measure  subtle  physiological  changes  in  the 
breast  that  are  indicative  of  disease  and  recovery  from  disease. 


III.  Progress  Summaries 


Ill.a  Clinical  Summary 

I  have  quantified  the  concentrations  of  blood,  water,  fat  and  scattering  in  the  breasts  of  over  35 
healthy  volunteers.  The  addition  of  fat  and  the  spectral  scattering  dependence  has  provided  some 
unique  insights  into  breast  physiology.  Current  data  have  demonstrated  : 
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•  NIR  absorption  and  scattering  in  the  pre-menopausal  breast  is  signifieantly  higher  than  in  the 
post-menopausal  breast  (Paper  #3,  Figures  2  &  3). 

•  Blood  and  water  decreases  as  age  increases  (Paper  #3,  Figures  4&  5). 

•  There  is  a  general  increase  in  lipids  as  age  increases  (Paper  #1 ,  Figure  4). 

•  Increased  scattering  in  the  pre-menopausal  breast  is  likely  due  to  higher  eollagen  density  and 
lower  amounts  of  fat  (the  eomponents  of  mammographic  density)  (Paper  #1,  Figure  6). 

•  Measurements  of  tissue  scattering  spectra  and  lipids  correlate  with  expected  breast 
eomposition  from  body-mass-index  (Paper  #1,  Figures  5  &6). 

•  Hemoglobin  saturation  (St02)  measurements  of  pre-menopausal  women  are  statistically 
lower  than  those  of  post-menopausal  women.  Post-menopausal  women  on  hormone 
replacement  therapy  (HRT)  show  St02  values  closer  to  the  pre-menopausal  breast  than  the 
post-menopausal  breast.  (The  St02  measures  of  the  arterial  and  venous  hemoglobin  saturation 
in  tissue)  (Paper  #3,  Table  2). 

•  Optical  assessment  of  tissue  physiology  is  sensitive  to  menstrual  cycle  (Paper  #2,  Table  1). 

All  of  the  findings  are  not  surprising,  but  until  now,  none  of  these  compositional  breast 
properties  of  the  tissue  could  be  measured  non-invasively  in  a  safe  manner.  These  measurements 
provided  a  control  to  which  we  could  compare  our  measurements  against  known  physiology.  At 
present,  they  represent  the  most  complete  set  of  data  available  on  the  NIR  optical  properties  of 
breast  tissue.  The  normal  subjects  provided  a  continuous  stream  of  patients  from  which  we  drew 
considerable  clinical  experience.  Detailed  studies  of  cancer  patients  did  not  commence  until 
early  2001  when  the  new  instrument  arrived  at  the  UCI  medical  center.  I  will  now  discuss  the 
results  of  several  exciting  measurements  performed  on  diseased  patients. 

One  published  measurement  clearly  demonstrates  our  ability  to  measure  important  physiological 
changes  brought  on  by  the  progression  of  disease  (Paper  #1,  Figures  11,12,  &13).  Measurements 
of  a  fibroadenoma  show  significant  increases  in  water  (Figure  11)  and  hemoglobin  (Figure  12). 
Traditional  optical  measures  of  tumors  have  relied  on  only  relative  measurements  (i.e.,  St02), 
whieh  was  slightly  lower  in  the  fibroadenoma  (Figure  13).  These  increases  in  water  and 
hemoglobin  have  also  been  seen  in  3  cancer  patients  (unpublished  data).  Interestingly,  the  St02 
was  significantly  lower  in  tumors  vs.  normal  compared  to  benign  lesions  vs.  normal 
(unpublished  results).  The  small  number  of  subjects  involved  does  not  provide  enough  statistical 
power  to  claim  this  will  always  be  the  case,  but  produces  a  physiological  basis  on  which  to 
proceed  for  continued  studies. 

Monitors  of  therapeutic  intervention  are  of  significant  interest  to  the  medical  community.  We 
have  measured  what  is  believed  to  be  the  first  optical  measurement  of  the  effects  of 
chemotherapy  in  a  human  subject  (to  be  published).  Our  results  demonstrate  : 
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•  For  a  responder,  tumor  water  and  hemoglobin  concentrations  decreased  significantly  over  the 
course  of  the  first  chemotherapy  cycle,  and  in  subsequent  cycles  as  well. 

•  For  a  responder,  normal  tissue  hemoglobin  concentrations  also  decreased  to  a  lesser  degree 
over  the  course  of  therapy.  Independent  measurements  of  hematocrit  (part  of  patient’s 
clinical  treatment,  and  not  our  trial)  also  decreased. 

•  For  a  responder,  fat  increased  and  tumor  oxygen  metabolism  (St02)  decreased  over  the 
course  of  therapy,  indicating  decreases  in  angiogenic  activity  of  the  tumor. 

•  The  final  lesion  prior  to  surgery  was  not  palpable  but  was  detectible  by  our  optical  method. 

•  Measurements  on  a  non-responder  showed  no  significant  changes  in  water  or  hemoglobin. 

These  preliminary  measurements  suggest  that  it  may  be  possible  to  monitor  the  effectiveness  of 

chemotherapy  using  a  non-invasive  technique.  Most  importantly,  it  is  clear  that  our  technique 

was  measuring  changes  in  tumor  angiogenesis  resulting  from  chemotherapy. 


IH.b  Instrumentation  Summary 

Two  technical  improvements  have  forwarded  the  achievement  of  the  initial  specific  aims.  We 
have  incorporated  a  steady-state  spectrometer  into  the  FDPM  instrument.  The  application  of  a 
simple  model  then  transforms  a  seven-wavelength  instrument  to  a  hundreds  of  wavelengths 
instrument.  The  spectrometer  added  previously  unimaginable  precision  in  determining  the 
physiological  makeup  of  tissues. 

I  have  also  completed  a  second  FDPM  instrument  that  resides  in  the  Breast  Center  at  the  Chao 
Comprehensive  Cancer  Center.  This  instrument  represents  significant  technical  advances, 
including  higher  spectral  content  and  faster  measurement  times.  This  new  system  was  used  in 
the  chemotherapy  measurements  listed  above  and  has  generated  physician  interest  in  the  project. 
The  success  of  this  instrument  led  us  to  complete  a  third  system. 


HI.  c  Protocol  Design/  Clinical  Exposure  Summary 

I  have  designed  and  managed  two  UCI  research  protocols  (#95-563  &  #99-2183).  I  have  added 
some  modifications  to  these  protocols  have  helped  us  gain  further  insight  into  the  optics  of  breast 
tissue.  No  difficulties  or  complications  have  occurred  during  the  course  of  these  trials. 
Exposure  to  this  aspect  of  clinical  research  has  taught  me  how  to  work  with  review  committees 
and  how  to  design  and  implement  clinical  research  protocols,  thus  improving  my  translational 
research  skills. 


Ill.d  Clinical  Collaboration 

The  close  connection  between  basic  and  clinical  research  was  one  of  the  primary  reasons  for 
coming  to  the  Beckman  Laser  Institute.  Interactions  with  clinicians  at  UCI  have  been  very 
finitfiil.  My  primary  physician  collaborators  have  been  Drs.  John  Butler  and  David  Hsiang,  both 
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surgical  oncologists.  Both  physicians  have  recruited  cancer  patients  and  provided  important 
medical  insight  into  our  analysis.  In  addition,  our  work  has  enabled  collaboration  with 
physicians  at  other  universities  (see  FV.d).  Their  interest  in  research  has  proven  invaluable  to  my 
development  as  a  breast  cancer  researcher. 


Ill.e  Student  Supervision 

I  have  also  benefited  from  working  with  Dr.  Tromberg’s  students.  I  have  been  partially 
responsible  for  training  and  guiding  them  in  their  studies  of  biomedical  optics.  Ryan  Lanning 
has  contributed  to  the  project  by  helping  construct  and  refine  instrumentation  used  in  our  breast 
studies,  and  has  now  graduated  and  gone  on  to  Harvard  Medical  School.  Natasha  Shah  has  been 
very  important  in  performing  clinical  measurements.  She  has  graduated  with  a  Master’s  in 
chemistry,  and  is  now  a  staff  researcher  in  Dr.  Tromberg’s  lab.  Dorota  Jakubowski,  an  MD/PhD 
student,  has  been  instrumental  as  a  programmer,  and  has  done  a  significant  amount  of  work  in 
the  chemotherapy  measurements  listed  above. 


IV  Conclusions 

Two  carefully-designed  experiments  demonstrated  that  the  FDPM  technique  measures 
functionally  relevant  physiological  changes  in  breast  tissue.  One  experiment  on  normal  women 
demonstrated  that  NIR-derived  physiology  correlates  nicely  with  known  histology.  The  other 
measurement  succeeded  in  monitoring  changes  in  tumor  angiogenesis.  Both  measurements  were 
performed  non-invasively.  These  results  now  form  a  foundation  upon  which  we  will  further 
refine  the  technique  and  push  focused  translational  research  applications  into  clinical  practice. 


V.  Appendix 

V.a  Reportable  Outcomes:  Peer-Reviewed  Publications 

[1]  A.  E.  Cerussi.  D.  Jakubowski,  N.  Shah,  F.  Bevilacqua,  R.  Lanning,  A.  J.  Berger,  D.  Hsiang,  J.  Butler,  R.  F. 
Holcombe,  and  B.  J.  Tromberg,  “Spectroscopy  enhances  the  information  content  of  optical 
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Diagnostic  methods  currently  in  use,  such  as  mammogra¬ 
phy,  magnetic  resonance  imaging,  and  ultrasound  (US), 
offer  excellent  anatomic  lesion-detection  capabilities  but 
generally  cannot  provide  quantitative  information  regard¬ 
ing  tissue  function  and  composition  (1).  Positron  emission 
tomography  shows  great  promise  in  evaluating  the  meta¬ 
bolic  demands  of  tissue  but  requires  exogenous  radionu¬ 
clides  and  is  insensitive  to  tissue  hemodynamics.  NIR 
optical  spectroscopy  is  intrinsically  sensitive  to  the  princi¬ 
pal  components  of  breast:  blood,  water,  and  adipose,  as 
well  as  epithelial  and  connective  tissues.  Results  of  pre¬ 
liminary  studies  suggest  that  the  fractional  contribution  of 
each  component  to  NIR  signals  depends  strongly  on  fac¬ 
tors  such  as  age,  menopausal  status,  and  the  progression 
of  disease.  Thus,  NIR  optical  spectroscopy  provides  an 
opportunity  for  revealing  physiologic  information  that 
cannot  otherwise  be  obtained  noninvasively. 
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As  optical  diagnostic  methods  become  more  clinically 
applicable,  there  is  a  growing  need  to  identify  and  quan¬ 
tify  sources  of  endogenous  contrast  in  breast  (2-1 1).  This 
article  presents  the  results  of  quantitative  NIR  spectro¬ 
scopic  studies  in  28  subjects  ranging  in  age  from  1 8  to  64 
years.  We  provide  evidence  that  NIR  tissue  spectroscopy 
is  inherently  sensitive  to  long-term  age  and  hormone-de¬ 
pendent  dynamics.  Similar  studies  have  been  initiated 
(12-14)  but  in  only  a  few  subjects  or  with  limited  sensi¬ 
tivity  to  physiologic  parameters.  These  data  are  expected 
to  play  an  important  role  in  characterizing  the  sensitivity 
of  NIR  optical  spectroscopy  and  imaging  for  applications 
such  as  tumor  detection,  therapeutic  monitoring,  and  risk 
assessment. 

NIR  photons  (600-1,000  nm)  are  nonionizing  and 
weakly  absorbed  by  tissue.  Many  investigators  have  used 
them  to  explore  brain,  breast,  and  mu.scle  tissue,  exploit¬ 
ing  their  enhanced  penetrance  relative  to  that  of  visible  or 
ultraviolet  radiation  (1).  In  breast  tissue,  the  principal  NIR 
absorbers  are  assumed  to  be  reduced  hemoglobin  (Hb-R), 
oxygenated  hemoglobin  (Hb-02),  water  (H2O),  and  lipids 
(15).  NIR  tissue  spectroscopy  offers  a  safe  way  to  quan¬ 
tify  these  components  and  to  view  unique  functional  in¬ 
formation  with  low-cost  “point-of-care”  devices.  Multiple 
light  scattering  complicates  quantitative  measurements  of 
light  absorption  in  the  NIR  range.  The  process  of  using 
frequency-domain  photon  migration  (FDPM)  techniques 
to  separate  the  effects  of  absorption  from  scattering  has 
been  described  in  detail  (16,17).  FDPM  enables  the  use  of 
quantitative  spectroscopic  analysis  tools  to  determine  the 
composition  and  structure  of  tissue  noninvasively  (18). 

The  application  of  FDPM  techniques  to  NIR  tissue 
spectroscopy  differs  substantially  from  eonventional  bulk 
NIR  methods,  such  as  diaphanography  (19).  FDPM  em¬ 
ploys  computational  models  that  allow  the  quantitative 
separation  of  the  effects  of  absorption  and  scattering  (20). 
Thus,  within  the  limits  of  applicability  of  the  light  trans¬ 
port  models  to  the  breast,  all  parameters  measured  in  this 
study  are  absolute  values. 


MATERIALS  AND  METHODS 


Patient  Selection 

All  volunteers  enrolled  in  this  study  competently  pro¬ 
vided  informed  consent  to  participate  in  one  of  two  trials 
(nos.  95-563  and  99-2183)  under  the  guidelines  of  an  Insti¬ 
tutional  Review  Board.  The  28  women  ranged  in  age  from 
18  to  64  years.  Fifteen  of  them  were  premenopausal  (aver¬ 
age  age,  28  years  ±  9),  and  seven  were  postmenopausal 


(average  age,  56  years  ±  2).  The  remaining  six  women  (av¬ 
erage  age,  56  years  ±  5)  were  taking  some  form  of  hormone 
replacement  therapy  (HRT);  three  of  the  six  classified  them¬ 
selves  as  perimenopausal.  None  of  the  women  had  any  known 
eancerous  lesions  or  other  known  forms  of  breast  disease. 

Measurement  Technique 

Each  subject  rested  in  a  supine  position  during  the  mea¬ 
surement.  The  instmment  probe,  slightly  larger  than  a  US 
probe  (Fig  1),  was  the  only  item  placed  in  contact  with  the 
subject.  It  consisted  of  a  sealed  APD  module  and  an  optical 
fiber  placed  22  mm  from  the  APD.  All  measurements  were 
performed  in  a  reflection-style  geometry.  The  probe  was 
placed  on  the  breast  with  minimal  pressure  and  only  the 
force  of  gravity;  no  compression  was  used.  In  this  configura¬ 
tion,  we  estimate  that  the  mean  sampling  depth  was  approxi¬ 
mately  1  cm  below  the  skin.  We  report  data  only  on  mea¬ 
surements  performed  in  the  center  of  the  left  upper  outer 
quadrant  (measurements  on  the  right  upper  quadrant  yielded 
similar  results).  Error  bars  in  the  figures  and  tables  represent 
the  standard  deviations  (SDs)  of  repeated  measurements. 

Instruments 

A  schematic  of  our  1-GHz,  portable  FDPM  device  is 
shown  in  Figure  1.  This  FDPM  instrument  has  been  de¬ 
scribed  in  detail  elsewhere  (21).  It  employs  multiple  di¬ 
ode  lasers  that  provide  visible  and  NIR  light  at  seven 
wavelengths  (672,  800,  806,  852,  896,  913,  and  978  nm). 
A  hand-held  probe  has  been  designed  to  house  an  APD 
that  records  the  modulated  diffuse  light  signals  after  prop¬ 
agation  through  the  tissue.  This  probe  has  a  plastic  attach¬ 
ment  on  the  casing  to  position  a  source  optical  fiber  a 
fixed  distance  from  the  APD.  A  100-p,m-diameter  graded 
index  optical  fiber  positioned  22  mm  from  the  APD  de¬ 
tector  was  used  to  deliver  the  diode  laser  outputs  to  the 
tissue  surface.  A  network  analyzer  (8753c;  Agilent,  Palo 
Alto,  Calif),  measured  the  phase  and  amplitude  of  the 
detected  electronic  signal  from  the  APD  over  a  range  of 
source  modulation  frequencies.  A  steady-state  current 
source  was  mixed  with  radiofrequency  power  provided  by 
the  network  analyzer  in  a  bias  network.  This  bias  network 
serially  distributed  power  to  each  laser  diode  and  pro¬ 
duced  intensity-modulated  light.  An  optical  switch  deliv¬ 
ered  light  serially  from  each  diode  to  the  tissue  via  the 
single  optical  fiber  described  above.  The  optical  power- 
launched  into  the  tissue  ranged  from  5  to  25  mW  for  each 
wavelength. 

A  sweep  over  all  seven  wavelengths  took  approximately 
35-60  seconds.  The  system  acquired  data  in  less  than  3  sec- 


212 


► 


Figure  1.  Schematic  drawing  of  FDPM  instrument,  hand-heid  probe,  and  measure¬ 
ment  map.  The  components  of  the  instrument  are  the  diode  iasers,  avaianche  photo¬ 
diode  (APD),  network  anaiyzer,  steady-state  current  source  (DC),  bias  network,  and 
opticai  switch.  Measurements  were  performed  with  a  hand-heid  probe  (ciose-up  at  bot¬ 
tom  of  Figure)  that  contained  the  APD  and  a  source  opticai  fiber.  The  measurement  ioca- 
tion  was  in  the  center  of  the  ieft  upper  outer  quadrant  of  the  breast,  as  indicated. 


onds  per  wavelength,  with  a  2-second  delay  between  each 
laser  diode  in  the  system  because  of  switching  consider¬ 
ations.  The  system  was  wheeled  into  a  medical  clinic  for 
each  measurement.  Instmmental  artifacts  were  removed  by 
means  of  calibration  on  a  tissue-simulating  phantom  with 
known  absorption  and  scattering  properties. 

Measured  Parameters  and  Data  Analysis 

The  amplitude  and  phase  of  a  NIR  diffusive  light  wave 
demodulates  and  retards,  respectively,  as  the  wave  propa¬ 


gates  through  multiple-scattering  media  such  as  tissue. 
The  real  and  imaginary  parts  of  this  diffuse  wave  were  fit 
simultaneously  to  a  light-diffusion  model  (the  Pi  approxi¬ 
mation  to  the  Boltzmann  transport  equation)  (22,23)  by 
minimizing  the  surface  with  a  Marquardt-Levenberg 
algorithm.  This  fit  determined  the  absolute  optical  absorp¬ 
tion  coefficient  ju-^  and  the  absolute  optical  reduced  scat¬ 
tering  coefficient  ju,'  at  each  wavelength.  Typical  jXa  and 
ju,^  uncertainties,  determined  from  the  distribution  of 
the  fits,  are  3%  to  1%  of  their  mean  values.  Precision 
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errors  in  actual  tissue  measurements  were  determined  by 
performing  multiple  network  analyzer  sweeps  and  multi¬ 
ple  probe  placements  on  a  given  location  and  were  found 
to  be  less  than  10%  for  ii,,  and  5%  for 

When  the  optical  properties  /a,'  and  ju,„  are  recovered 
for  the  seven  wavelengths,  the  spectral  dependence  of  the 
absorption  may  be  combined  with  known  values  of  mo¬ 
lecular  extinction  coefficients  to  calculate  physiologically 
relevant  parameters.  We  assume  that  bulk  breast  tissue  is 
composed  principally  of  four  NIR  absorbers:  Hb-R, 

Hb-02,  H2O,  and  lipids.  The  concentrations  of  these  ab¬ 
sorbers  (c)  were  quantified  by  solving  the  equation  = 
Ec,  where  £  is  a  7  X  4  matrix  that  contains  the  molar 
extinction  coefficients  of  the  four  chromophores  at  the 
seven  source  wavelengths  (24-26). 

For  each  measurement,  we  report  four  hemoglobin  pa¬ 
rameters:  [Hb-R],  [Hb-Oi],  the  total  hemoglobin  concen¬ 
tration  (THC  =  [Hb-R]  -I-  [Hb-Oi]),  and  the  tissue  hemo¬ 
globin  saturation  (SiOi  =  [Hb-O^i/THC  X  100%),  where 
the  characters  in  square  brackets  denote  concentration  in 
moles  per  liter.  Values  for  water  and  lipid  content  are 
reported  as  percentages.  The  water  percentage  is  the  con¬ 
centration  of  measured  tissue  water  divided  by  the  con¬ 
centration  of  pure  water  (55.6  mol/L).  The  lipid  percent¬ 
age  indicates  the  kilograms  per  liter  of  lipid  measured 
relative  to  an  assumed  “pure”  lipid  density  of  0.9  kg/L. 
Thus,  the  water  and  lipid  percentages  we  report  are  rela¬ 
tive  figures  of  merit  compared  with  pure  solutions  of  the 
substance  and  are  not  strict  volume  or  mass  fractions. 

The  scattering  properties  of  the  tissue  also  yield  important 
physiologic  information.  NIR  scattering  in  tissue  has  the 
following  dependence:  where  A  is  a  constant,  A 

is  the  wavelength  (in  nanometers),  and  SP  is  the  scatter 
power.  Scatter  power  is  related  to  the  scattering  center  size 
(d)  compared  with  the  optical  wavelength.  As  an  example, 
scatter  power  is  4  in  the  case  of  Rayleigh  scattering  (d  ^  A) 
and  is  approximately  1  for  large  Mie-like  scatterers  {d  ~  A). 


RESULTS 


Sensitivity:  Pre-  versus  Postmenopausal  Subjects 

Both  Figures  2  and  3  present  a  typical  series  of  mea¬ 
surements  of  seven  absorption  and  scattering  coefficients, 
respectively.  The  points  represent  an  average  of  several 
measurements  in  the  center  of  the  left  upper  outer  breast 
quadrants  in  two  subjects,  a  32-year-old  premenopausal 
woman  and  a  54-year-old  postmenopausal  woman.  Error 
bars  indicate  the  SDs  of  repeated  measurements.  There 
are  vast  absorption  and  scattering  differences  between 


Wavelength  (nm) 

Figure  2.  Measured  absorption  spectra  for  a  32-year-old  pre¬ 
menopausal  subject  (■,  solid  line)  and  a  54-year-old  postmeno¬ 
pausal  subject  (A,  dashed  line).  Points  represent  the  average  of 
several  measurements,  and  lines  represent  a  least-squares  fit  (ex¬ 
trapolated  to  all  wavelengths)  based  on  the  assumption  that 
breast  absorption  is  due  only  to  Hb-R,  Hb-Oz,  H2O,  and  lipids. 
Error  bars  show  the  SDs  of  repeated  measurements.  Concentra¬ 
tions  of  chromophores  from  this  measurement  are  listed  in  Table  1 . 


Figure  3.  Measured  scattering  spectra  for  the  same  subjects  as 
in  Figure  2.  The  lines  represent  a  fit  to  =  AA”®'’,  where  fXs  is 
the  reduced  scattering  coefficient  (per  millimeter),  A  is  a  constant, 
A  is  the  wavelength  (in  nanometers)  and  SP  is  the  scatter  power. 
Results  of  the  fit  from  this  measurement  are  listed  in  Table  1.  Er¬ 
ror  bars  show  the  SDs  of  repeated  measurements. 

pre-  and  postmenopausal  breast  tissue.  The  solid  lines  in 
Figure  2  represent  a  weighted  least-squares  fit  of  the 
seven  absorption  coefficients,  with  published  extinction 
coefficients  used  for  Hb-R,  Hb-02,  HiO,  and  lipids  (24- 
26).  Lines  between  the  measured  points  have  been  inter¬ 
polated.  The  solid  lines  of  Figure  3  represent  a  fit  of  the 
scattering  coefficients  to  the  equation  =  AA^^’’. 
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NIR  OPTICAL  CONTRAST  IN  BREAST 


Table  1 

Measured  Physiologic  Properties  in  Breasts  of  a  Premenopausal  and  a  Postmenopausal  Subject 


Subject 

[Hb-R]  (/j,mol/L) 

[Hb-02]  (/J.mol/L) 

S,02  {%) 

THC  (/xmol/L) 

[H2O]*  (%) 

Lipidt  (%) 

Scatter  Power 

Premenopausal 

Postmenopausal 

12.6  ±  0.7 

2.4  ±  0.5 

27.9  ±  2.7 

12.2  ±  1.6 

68.9  ±  1.3 

83.7  ±  2.5 

40.4  ±  2.7 

14.4  ±1.9 

45.1  ±  3.6 

10.3  ±  0.8 

29.8  ±  5.6 

67.5  ±  3.7 

0.864  ±  0.068 

0.555  ±  0.036 

Note. — Measurements  are  given  as  individual  means  ±  SDs  for  the  two  subjects  represented  in  Figures  2  and  3. 
*  Measured  concentration  relative  to  pure  water, 
t  Measured  mass  density  relative  to  pure  lipid. 


Figure  4.  Plot  of  the  THC  versus  subject  age.  Premenopausal 
subjects  (<50  years,  n  =  15)  show  considerable  variation  with 
age  because  of  intersubject  variations,  such  as  menstrual  cycle 
variations  and  differences  in  overall  hormone  production.  Post¬ 
menopausal  subjects  (>50  years,  n  =  1 3)  appear  to  show  a  de¬ 
crease  in  THC.  Error  bars  show  the  SDs  of  repeated  measure¬ 
ments. 

These  optical  spectra  provide  insight  into  breast  tissue 
composition.  There  are  higher  concentrations  of  hemoglo¬ 
bin  (ie,  both  Hb-R  and  Hb-02)  in  the  premenopausal  sub¬ 
ject,  as  evidenced  by  the  overall  higher  absorption  in  the 
670-850-nm  range.  There  is  also  more  water  relative  to 
lipids  in  the  premenopausal  subject,  as  revealed  by  the 
large  water  absorption  peak  at  980  nm.  Recovery  of  the 
absorption  spectra  allows  calculations  of  the  tissue  con¬ 
centrations  of  Hb-R,  Hb-02,  H2O,  and  lipids.  In  addition, 
the  light  scattering  intensity  is  substantially  lower  in  the 
postmenopausal  breast.  Table  1  summarizes  the  fitted 
physiologic  properties  for  the  two  subjects  represented  in 
Figures  2  and  3.  ^ 

Quantities  versus  Age 

The  next  three  figures  present  FDPM-measured  quanti¬ 
ties  as  a  function  of  subject  age.  Figure  4  shows  the 


THC,  Figure  5  shows  the  tissue  water  percentage  (refer¬ 
enced  to  the  concentration  of  pure  water),  and  Figure  6 
shows  the  scatter  power.  Premenopausal  subjects  (ie,  aged 
younger  than  50  years)  display  a  variety  of  values  in  all 
three  plots.  This  spread  is  the  result  of  intersubject  varia¬ 
tions,  including,  but  not  limited  to,  menstrual  cycle  differ¬ 
ences  and  gynecologic  age.  The  THC,  and  to  a  lesser  ex¬ 
tent  the  water,  seemed  to  increase  in  premenopausal  sub¬ 
jects  (aged  18-39  years),  perhaps  reaching  a  peak  value 
near  age  30  years.  Between  the  ages  of  40  and  49  years, 
the  THC  appears  to  level  off,  while  water  and  scatter 
power  seem  to  decline.  After  age  50  years  (predominantly 
postmenopausal  subjects),  there  is  a  general  decrease  with 
age  in  THC,  water,  and  scatter  power 

The  late  decrease  in  THC  and  water  correlates  well 
with  previous  histologic  studies  showing  both  atrophy  in 
well-vascularized  lobular  tissue  and  an  increase  in  the 
fat-to-collagen  ratio  after  menopause  (15).  This  decrease 
is  consistent  with  compositional  analysis  data  showing 
lower  blood  and  water  content  for  fat  versus  glandular 
tissue  (27).  Scatter  power  also  declines  with  age  after  age 
50  years.  Compared  with  fat,  collagen  and  glandular  tis¬ 
sue  scatter  light  with  a  higher  intensity  and  a  steeper 
spectral  dependence  (28).  Thus,  smaller  scatter  powers  are 
expected  in  fatty  tissue.  Exact  knowledge  of  the  weighted 
biologic  contributions  to  optical  tissue  scattering  is  a  mat¬ 
ter  of  debate. 

Correlation  between  Absorption  and  Scattering 
Information 

Figure  7  displays  a  correlation  plot  of  the  measured 
water  and  lipid  percentages  versus  scatter  power.  The 
lipid  percentage  is  the  mass  density  of  lipids  in  the  tissue 
referenced  to  the  mass  density  of  soybean  oil.  Note  that 
neither  the  water  nor  the  lipid  percentage  represents  a 
strict  mass  or  volume  fraction  for  the  tissue  (ie,  not  all 
components  are  normalized  to  100%).  Lines  represent  fits 
of  water  (solid  line)  and  lipids  (dashed  line)  as  functions 
of  scatter  power. 


215 


CERUSSI  ET  AL 


Academic  Radiology,  Vol  8,  No  3,  March  2001 


Figure  5.  Plot  of  the  measured  tissue  water  percentage  versus 
subject  age.  Premenopausal  subjects  show  considerable  variation 
over  age  because  of  intersubject  variations  such  as  menstrual 
cycle  variations  and  differences  in  overall  hormone  production. 
Postmenopausal  and  perimenopausal  subjects  appear  to  show  a 
decrease  in  water  content.  The  percentages  are  water  concentra¬ 
tions  normalized  to  the  concentration  of  pure  water.  Error  bars 
show  the  SDs  of  repeated  measurements. 


Tissue  scattering  correlates  well  with  the  parenchymal 
composition  measured  via  absorption.  As  scatter  power 
increases,  the  lipid  percentage  decreases  (R-  =  -0.841, 
P  <  .001)  and  the  water  percentage  increases  (R-  = 
-1-0.850,  P  <  .001).  The  absorption  spectra  of  water 
and  lipids  overlap,  which  can  lead  to  errors  in  measuring 
these  chromophores.  Figure  7,  however,  clearly  shows 
that  both  scattering  and  absorption  provide  complemen¬ 
tary  information  regarding  breast  composition,  indicating 
that  such  errors  do  not  substantially  distort  our  findings. 

Tissue  Metabolism  and  Menopause 

The  S,02  is  the  percentage  fraction  of  arterial  and  ve¬ 
nous  Hb-Oi  that  contributes  to  the  tissue  THC,  and  it 
should  not  be  confused  with  arterial  saturation,  which  is 
measured  with  pulse  oximetry.  Table  2  presents  the  mea¬ 
sured  S,02  in  all  subjects  with  a  body-mass  index  less 
than  40.  The  average  S,02  was  approximately  76%  for 
both  the  premenopausal  subjects  and  those  receiving 
HRT;  their  means  were  not  significantly  different.  The 
S,02  values  for  the  postmenopausal  women  not  receiving 
HRT,  however,  were  significantly  higher,  averaging 
81.9%  ±  8.8  (P  =  .062). 

Decreases  in  S,02  are  associated  with  increased  metab¬ 
olism  (ie,  extraction  of  O2  from  the  blood).  Tumors,  for 


Figure  6.  Plot  of  the  measured  optical  scatter  power  (in  arbitrary 
units)  versus  subject  age.  Optical  scatter  power  refers  to  the 
magnitude  of  the  slope  of  the  scattering  versus  wavelength  curve. 
Premenopausal  subjects  show  considerably  higher  scatter  power 
than  postmenopausal  subjects.  This  general  pattern  is  the  same 
for  mammographic  density.  Error  bars  show  the  SDs  of  repeated 
measurements. 


Scatter  Power 


Figure  7.  Correlation  plots  for  water  percentage  (□,  solid  line) 
and  lipid  percentage  (A,  dashed  line)  versus  scatter  power.  Per¬ 
centages  are  referenced  to  the  concentration  of  pure  water  and 
the  mass  density  of  pure  fat.  Scatter  power  correlates  with  water 
(R2  =  +0.850,  P  <  .001)  and  correlates  negatively  with  lipids 
(P2  =  -0.841,  P  <  .001),  Error  bars  show  the  SDs  of  repeated 
measurements. 

example,  have  been  .shown  in  vivo  to  have  lower  SjOt 
values  than  nontumor  tissue  (6,18,29).  The  postmeno¬ 
pausal  increase  in  S,02  relative  to  the  premenopausal 
women  and  those  receiving  HRT  is  plausible,  because 
general  tissue  metabolism  in  the  breast  should  decrease 
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NIR  OPTICAL  CONTRAST  IN  BREAST 


Table  2 

Comparison  of  St02  Measurements 


Subjects 

StOa  (%) 

Premenopausal  (n  =  15) 

76.3  ±  6.4 

Receiving  HRT  {n  =  6) 

76.9  ±  9.2 

Postmenopausai  (n  =  6) 

81.9  ±  8.8 

Note. — Measurements  are  given  as  means  ±  SDs  and  are  pro¬ 
vided  for  aii  subjects  with  a  body-mass  index  of  iess  than  40. 


after  menopause  as  the  tissue  begins  to  atrophy.  The  simi¬ 
larity  between  8(02  values  in  premenopausal  and  HRT 
women  suggests  that  their  metabolisms  are  similar  be¬ 
cause  of  the  effects  of  both  endogenous  and  exogenous 
hormones. 


DISCUSSION 


Although  it  is  difficult  to  validate  noninvasive  in  vivo 
measurements  directly,  our  initial  results  indicate  that  the 
sensitivity  of  NIR  spectroscopy  is  a  reasonable  reflection 
of  long-term  hormone-controlled  breast  remodeling.  Addi¬ 
tional  work  in  the  field  has  demonstrated  that  short-term 
menstrual  cycle  changes  are  detectible  with  a  similar  NIR 
technique  (30).  Our  results  indirectly  validate  the  general 
accuracy  of  NIR  breast  spectroscopy.  This  quantitative 
physiologic  information  is  not  obtainable  with  any  other 
noninvasive  radiologic  method. 

An  important  feature  of  NIR  methods  is  the  ability  to 
characterize  quantitatively  the  breast  tissue  of  women  re¬ 
gardless  of  age,  hormonal  status,  or  mammographic  den¬ 
sity.  Increased  mammographic  density  contributes  to  a 
22%  false-negative  rate  and  a  high  false-positive  rate 
(56.2%  cumulative  risk  after  10  examinations)  in  women 
less  than  50  years  of  age  (31,32).  A  recent  study  found 
that  routine  initial  mammography  was  not  clinically  ad¬ 
vantageous  for  women  less  than  35  years  of  age  (33). 
Furthermore,  the  use  of  HRT  in  postmenopausal  women 
is  known  to  increase  mammographic  density  (34)  and  has 
been  recently  shown  to  impede  the  efficacy  of  mammo¬ 
graphic  screening  (35,36).  The  sensitivity  of  NIR  spec¬ 
troscopy  to  known  biologic  processes  suggests  that  optical 
methods  may  reveal  important  new  information  comple¬ 
menting  that  provided  by  conventional  diagnostic  tech¬ 
niques,  particularly  in  the  case  of  radiographically  dense 
breast  tissue. 

In  general,  premenopausal  breast  tissue  is  more  opti¬ 
cally  attenuating  than  postmenopausal  breast  tissue,  in 


both  absorption  and  scattering.  It  is  important  to  note  that 
this  trend,  particularly  in  scattering,  is  also  true  for  mam¬ 
mographic  density.  Radiographically  dense  breast  tissue  is 
due  to  differing  amounts  of  fat,  collagen,  epithelium,  and 
water  (37).  Although  x  rays  are  inherently  sensitive  to 
dense  electron  structures  in  tissue,  such  as  microcalcifica¬ 
tions,  NIR  photons  are  naturally  sensitive  to  tissue  con¬ 
stituents.  Thus,  NIR  spectroscopy  has  a  unique  potential 
for  quantifying  the  elements  of  breast  tissue  that  contrib¬ 
ute  to  mammographic  density.  This  observation  may  be 
important  in  screening,  since  it  may  identify  breast  tissue 
at  physiologic  risk  for  malignant  transformation,  and  the 
technique  can  be  performed  easily  in  all  women. 

In  addition  to  scattering  contrast,  we  report  age-depen- 
dent  variations  in  THC,  H2O,  and  lipids,  as  well  as 
changes  in  8,02  for  pre-  and  postmenopausal  subjects. 
Knowledge  of  the  “normal”  values  of  these  chromophores 
will  help  in  evaluating  the  usefulness  of  optical  methods 
for  detecting  and  characterizing  lesions  in  the  breast.  8ev- 
eral  investigators  have  reported  a  two-  to  fourfold  THC 
contrast  between  normal  and  tumor  structures.  In  vivo 
tumor  8,02  values  are  also  typically  lower  than  those  of 
normal  tissue  (6,18,29).  Thus,  detailed  studies  of  normal 
tissue  are  essential  for  determining  the  sensitivity  required 
of  optical  instrumentation  for  detecting  lesions  in  women 
of  varying  age  and  hormonal  status.  Furthermore,  as  base¬ 
line  levels  are  characterized,  data  on  an  individual’s  ab¬ 
sorption  and  scattering  variations  could  provide  important 
insight  into  disease  appearance  and  progression.  When 
applied  to  patients  receiving  chemotherapy  and/or  HRT, 
this  information  could  also  be  used  to  generate  feedback 
that  would  permit  customized  treatment  planning  based 
on  individual  physiologic  response. 
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Near  infrared  diffuse  optical  spectroscopy  and  diffuse  optical 
imaging  are  promising  methods  that  eventually  may  enhance  or 
replace  existing  technologies  for  breast  cancer  screening  and 
diagnosis.  These  techniques  are  based  on  highly  sensitive,  quan¬ 
titative  measurements  of  optical  and  functional  contrast  between 
healthy  and  diseased  tissue.  In  this  study,  we  examine  whether 
changes  in  breast  physiology  caused  by  exogenous  hormones, 
aging,  and  fluctuations  during  the  menstrual  cycle  result  in  signif¬ 
icant  alterations  in  breast  tissue  optical  contrast.  A  noninvasive 
quantitative  diffuse  optical  spectroscopy  technique,  frequency- 
domain  photon  migration,  was  used.  Measurements  were  per¬ 
formed  on  14  volunteer  subjects  by  using  a  hand-held  probe. 
Intrinsic  tissue  absorption  and  reduced  scattering  parameters  were 
calculated  from  frequency-domain  photon  migration  data.  Wave¬ 
length-dependent  absorption  (at  674,  803,  849,  and  956  nm)  was 
used  to  determine  tissue  concentration  of  oxyhemoglobin,  deoxy¬ 
hemoglobin,  total  hemoglobin,  tissue  hemoglobin  oxygen  satura¬ 
tion,  and  bulk  water  content.  Results  show  significant  and  dra¬ 
matic  differences  in  optical  properties  between  menopausal  states. 
Average  premenopausal  intrinsic  tissue  absorption  and  reduced 
scattering  values  at  each  wavelength  are  2.5-  to  3-fold  higher  and 
16-28%  greater,  respectively,  than  absorption  and  scattering  for 
postmenopausal  subjects.  Absorption  and  scattering  properties  for 
women  using  hormone  replacement  therapy  are  intermediate 
between  premenopausal  and  postmenopausal  populations.  Phys¬ 
iological  properties  show  differences  in  mean  total  hemoglobin 
(7.0  fiM,  11.8  /xM,  and  19.2  p.M)  and  water  concentration  relative 
to  pure  water  (10.9%,  15.3%,  and  27.3%)  for  postmenopausal, 
hormone  replacement  therapy,  and  premenopausal  subjects,  re¬ 
spectively.  Because  of  their  unique,  quantitative  information  con¬ 
tent,  diffuse  optical  methods  may  play  an  important  role  in  breast 
diagnostics  and  improving  our  understanding  of  breast  disease. 

X-ray  mammography  is  widely  used  for  screening  breast 
cancer,  the  most  common  form  of  cancer  in  women.  Because 
of  uncertainties  associated  with  radiographic  density,  mammog¬ 
raphy  has  up  to  a  22%  false-negative  rate  as  well  as  a  high 
false-positive  rate  (56.2%  cumulative  risk  after  10  exams)  in 
women  under  50  years  of  age  (1,  2).  A  recent  study  found  that 
routine  initial  mammography  was  not  clinically  advantageous  for 
women  under  35  years  of  age  (3).  Furthermore,  the  use  of 
hormone  replacement  therapy  (HRT)  in  postmenopausal 
women  is  known  to  increase  mammographic  density  (4)  and  has 
been  shown  recently  to  impede  the  efficacy  of  mammographic 
screening  (5,  6).  Techniques  such  as  MRI  and  ultrasound  (US) 
are  used  only  as  secondary  procedures  because  of  factors  such  as 
high  cost  and  poor  specificity  (MRI)  or  low  sensitivity  (US). 

Currently,  invasive  procedures  such  as  fine-needle  aspiration 
or  surgical  biopsy  are  implemented  to  provide  a  definitive 
diagnosis.  Given  the  suboptimal  performance  of  x-ray  mammog¬ 
raphy  in  premenopausal  and  perimenopausal  women,  the  ma¬ 
jority  of  invasive  follow-up  procedures  are  performed  on  normal 
or  benign  tissue  that  present  no  malignant  disease  (7).  As  a 
result,  the  use  of  noninvasive,  near-infrared  (NIR)  optical 


methods  as  a  supplement  to  present  techniques  for  diagnosing 
and  detecting  breast  cancer  has  generated  considerable  interest. 

Optical  methods  are  advantageous,  because  they  are  nonin¬ 
vasive,  quantitative,  and  relatively  inexpensive;  do  not  require 
compression;  and  pose  no  risk  of  ionizing  radiation.  A  promising 
NIR  optical  technique  currently  under  development  is  frequency 
domain  photon  migration  (FDPM).  FDPM  methods  have  de¬ 
tected  successfully  the  presence  of  small  palpable  breast  lesions 
in  vivo  in  women  with  previously  diagnosed  breast  abnormalities 
(8-10). 

FDPM  employs  intensity-modulated  NIR  light  to  characterize 
tissue  quantitatively  in  terms  of  its  optical  parameters,  i.e.,  the 
reduced-scattering  coefficient  (p,s')  and  the  absorption  coeffi¬ 
cient  (/Lia).  The  concentration  of  significant  NIR  absorbers 
[deoxyhemoglobin  (Hb),  oxyhemoglobin  (Hb02),  water,  and 
adipose]  can  be  calculated  by  using  measured  jUa  values  (11, 12). 
Multiple  scattering  of  light  in  breast  tissue  occurs  as  a  conse¬ 
quence  of  spatial  variations  in  refractive  index,  which  are  influ¬ 
enced  by  cellular  and  extracellular  matrix  density  (13,  14).  Both 
jHa  and  p.s'  provide  an  understanding  of  changes  in  tissue 
cellularity,  metabolic  activity,  physiology,  and  host  response 
to  cancer.  Detection  of  lesions  is  based  on  the  functional  contrast 
between  normal  and  diseased  tissue  in  the  same  patient. 
However,  the  physiology  of  healthy  breast  tissue  is  complex, 
influenced  by  multiple  internal  and  external  factors  such  as 
menstrual-cycle  phase,  menopausal  status,  exogenous  hor¬ 
mones,  lactation,  and  pregnancy.  Consequently,  to  establish  a 
better  basis  for  optical  detection  and  diagnosis  based  on  differ¬ 
ential  functional  contrast,  the  optical  properties  of  normal  breast 
tissue  must  be  carefully  examined  and  characterized. 

There  has  been  increasing  attention  paid  to  the  study  of 
normal  breast  tissue  optical  properties  and  changes  that  occur 
with  age,  menopausal  status,  and  HRT  (15-17).  During  the 
reproductive  years,  the  breast  is  comprised  mostly  of  glandular 
tissue.  With  menopause,  there  is  a  progressive  atrophy  of 
glandular  tissue;  mitotic  activity  slows  down;  and  vascular  re¬ 
quirements  decrease.  Some  of  these  changes  have  been  moni¬ 
tored  qualitatively  by  mammography  and  MRI  (3,  18-20). 
Involution  of  breast  glandular  tissue  can  be  impeded  by  the  use 
of  exogenous  hormones,  such  as  estrogen  and  progesterone  (21). 
Presently,  up  to  20%  of  U.S.  women  use  HRT  to  alleviate 
symptoms  of  menopause.  However,  long-term  use  has  been 
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linked  with  increased  risk  for  breast  cancer  (22).  MRl  studies 
further  sugge.st  that  fluctuations  in  endogenous  hormones  may 
alter  menstrual-  cycle  patterns  and  impact  breast  cancer  risk 
(23-26).  In  addition,  the  phase  of  the  menstrual  cycle  has  been 
shown  to  influence  mammographic  accuracy  and  survival  rates 
after  tumor  resection  (27,  28). 

In  this  article,  we  provide  initial  results  examining  differences 
in  both  optical  and  physiological  properties  of  healthy  breast 
tissue  of  14  women.  Our  results  show  that  the  population  studied 
can  be  well  differentiated  on  the  basis  of  quantitative  NIR 
spectroscopic  measurements.  Optical  contrast  is  based  on  phys¬ 
iological  changes  accompanying  menstrual  cycle  variations,  ex¬ 
ogenous  hormone  levels,  and  menopausal  status.  We  believe  this 
study  is  particularly  important,  because  standard-of-care  .screen¬ 
ing  with  x-ray  mammography  has  diminished  efficacy  in  two  of 
the  populations  we  examine,  i.e.,  premenopausal  and  HRT 
subjects.  The  sensitivity  of  our  technique  to  known  biological 
processes  suggests  our  methods  may  provide  important  infor¬ 
mation  complementary  to  conventional  diagnostic  techniques, 
particularly  in  the  case  of  radiographically  dense  breast  tis.sue. 
Ultimately,  this  technique  may  enhance  our  understanding  of 
pathophysiological  changes  that  accompany  malignant  transfor¬ 
mation  and  provide  insight  into  processes  associated  with  in¬ 
creased  disease  risk. 

Materials  and  Methods 

FORM  Instrument.  A  portable,  multiwavelength,  high-bandwidth 
FDPM  instrument  has  been  designed  and  optimized  for  clinical 
optical  property  studies  (Fig.  1,  quadrant  a).  The  instrument 
employs  multiple  diode  lasers  (box  1 )  to  provide  visible  and  NIR 
light  at  six  wavelengths  (674,  803,  849,  894,  947,  and  956  or  980 
nm).  The  FDPM  instrument  is  described  in  detail  by  Tromberg 
et  al.  (10)  and  Pham  et  al.  (29).  A  hand-held  probe  has  been 
designed  to  house  an  avalanche  photodiode  (A  PD)  that  records 
the  diffuse  light  signals  after  propagation  through  the  tissue  (box 
2).  This  probe  has  multiple  grooves  on  the  casing  to  position 
source  optical  fibers  a  fixed  distance  from  the  APD.  A  100-/Lim- 
diameter  graded  index  optical  fiber  positioned  25  mm  from  the 
APD  detector  delivers  the  diode  laser  output  to  the  tissue 
surface.  The  network  analyzer  (box  3)  measures  the  phase  and 
amplitude  of  the  electronic  signal.  A  dc  source  (box  4)  is  mixed 
with  rf  power  provided  by  the  network  analyzer  in  a  bias  network 
(box  5),  which  distributes  power  to  each  laser  diode  and  produces 
amplitude  modulated  light.  An  optical  switch  (box  6)  delivers 
light  serially  from  each  diode  to  the  tissue.  The  optical  power 
launched  into  the  subject  ranged  from  10  to  20  MW  for  each 
wavelength. 

Measurement  Technique.  Subjects  were  measured  in  the  supine 
position.  Measurements  were  made  in  reflectance  geometry  by 
placing  the  hand-held  probe  on  the  skin  surface  with  light 
pressure.  The  phase  shift  and  amplitude  of  the  photon  density 
wave  was  measured  at  each  wavelength  for  201  modulation 
frequencies  between  50  MHz  and  1,000  MHz  at  a  source- 
detector  separation  of  25  mm.  The  range  of  modulation  fre¬ 
quencies  was  swept  repetitively  such  that  each  amplitude  and 
phase  value  represents  up  to  12  measurements.  Four  to  six 
wavelengths  were  used  for  the  measurements. 

Subjects.  A  total  of  14  subjects  was  studied:  Pre  subjects  (Pre, 
premenopausal)  1-6  were  premenopausal  subjects  aged  18  to  33 
years.  Pre  subject  4  (29  years  old)  was  measured  during  the 
proliferative  phase  of  her  menstrual  cycle;  Pre  subject  1  (27  years 
old)  was  measured  around  ovulation;  and  Pre  subjects  2,  3,  and 
5  (aged  22,  26,  and  32  years,  respectively)  were  measured  during 
the  secretory  phase.  Pre  subject  6  (18  years  old)  was  studied  at 
two  different  points  in  her  menstrual  cycle.  One  measurement 
was  performed  near  ovulation  (day  14)  of  the  subject’s  regular 


Fig.  1.  Schematic  drawing  of  FDPM  instrument,  hand-held  probe,  and 
measurement  map  of  healthy  subjects.  The  components  of  the  instrument  are: 
diode  lasers  (box  1),  avalanche  photodiode  (box  2),  network  analyzer  (box  3), 
DC  current  source  (box  4),  bias  T  (box  5),  and  optical  switch  (box  6).  See  text  for 
detailed  description.  The  breast  is  divided  into  four  quadrants:  upper  outer 
(a),  upper  inner  (b),  lower  outer  (c),  and  lower  inner  (d).  FDPM  measurements 
are  made  in  each  quadrant,  on  the  areolar  border  (e),  and  on  the  glandular  tail 
(f)  that  extends  into  the  axilla. 

28-day  cycle,  and  another  measurement  was  performed  during 
the  luteal  phase  (day  25)  of  the  menstrual  cycle.  Post  subjects 
(Post,  postmenopausal)  1-3  were  postmenopausal  women  aged 
57  to  67  years.  The  remaining  five  subjects  were  postmenopausal 
women  51-60  years  of  age  who  were  using  HRT.  HRT  subjects 
1  and  5  were  on  combination  HRT  (progesterone  and  estrogen), 
and  HRT  subjects  2-4  were  taking  estrogen-only  HRT.  All  HRT 
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subjects  have  been  taking  exogenous  hormone  for  at  least  2 
years.  Experiments  were  conducted  in  adherence  to  University 
of  California,  Irvine,  Institutional  Review  Board-approved  pro¬ 
tocols  95-563  and  99-2183.  The  subjects  were  healthy  with  no 
known  breast  diseases.  After  providing  informed  consent,  sub¬ 
jects  filled  out  a  brief  questionnaire  that  surveyed  pertinent 
medical  history. 

When  using  four  wavelengths,  the  time  required  for  a  single 
measurement  was  approximately  90  seconds.  A  total  of  12  sites, 
six  per  breast,  for  each  subject  was  measured  during  subject 
visits.  The  measurement  map  is  shown  in  Fig.  1.  Measurements 
were  made  in  each  of  the  four  quadrants,  approximately  halfway 
from  the  center  of  the  breast  to  the  edge  of  the  breast,  laterally 
(a-d).  Additional  measurements  were  made  on  the  superior 
areolar  border  (e)  and  the  axillary  tail  (f).  Calibration  measure¬ 
ments  were  made  approximately  every  10  minutes  on  a  tissue 
phantom  of  known  optical  properties. 

Model.  The  Pi  approximation  to  the  Boltzmann  transport  equa¬ 
tion  (30)  was  used  to  analyze  the  frequency-domain  data  and  is 
valid  when  probing  homogeneous  turbid  media  at  source- 
detector  separations  greater  than  1  cm  and  modulation  frequen¬ 
cies  less  than  a  few  GHz.  The  model  uses  an  extrapolated 
boundary  condition  for  a  semiinfinite  geometry  (31).  The  solu¬ 
tion  to  the  P\  approximation  relates  the  frequency-dependent 
phase  and  amplitude  to  the  /Xa  and  /j-s'  coefficients  (32).  To 
determine  the  optical  properties  from  a  given  set  of  frequency- 
dependent  data,  a  Marquardt-Levenburg  minimization  algo¬ 
rithm  was  used  to  fit  simultaneously  the  amplitude  A{u>)  and 
phase  4>(c<j)  by  minimizing  the  difference  between  the  measured 
values  and  those  predicted  by  the  Pi  approximation. 

Physiological  properties  were  calculated  from  the  determined 
jXa  values  assuming  principal  chromophores  in  breast  tissue  for 
the  NIR  wavelengths  used  are  Hb,  Hb02,  and  water.  The 
concentrations  in  tissue  of  the  three  components  were  calculated 
by  using  p,a  at  four  wavelengths:  674, 803, 849,  and  956  or  980  nm. 
The  contribution  of  each  component  to  the  total  /Xa  at  a  given 
wavelength  is  represented  by  the  equation  2.303(ei  cj)  =  /Xai, 
where  ei  is  the  extinction  coefficient  in  units  of  cm^/mol  of  a 
particular  chromophore  at  a  given  wavelength  (A)  and  Ci  is  the 
concentration  of  the  chromophore.  The  total  absorption  at  a 
given  wavelength  from  the  tissue  chromophores  is 

s\Hb)[Hb]  +  e\Hboj[Hb02]  +  e\H,0)[H20]  =  [1] 

where  the  brackets  ([  ])  indicate  concentration.  The  matrix 
representation  of  the  four  equations  is  (11) 


Frequency  [MHz] 


Fig.  2.  FDPM  measurements  of  phase  lag  (a)  and  amplitude  vs.  modulation 
frequency  (b)  for  Pre  1  and  Post  1 .  Source-detector  separation,  2.5  cm;  wave¬ 
length,  849  nm.  Solid  lines  represent  best  diffusion-model  function  fits  to 
phase  and  amplitude  data  for  each  subject.  Error  bars  are  on  the  order  of  the 
marker  size. 


To  extract  4>niedium  from  <I>meas  at  a  single  source-detector 
separation,  '^instrument  must  be  explicitly  evaluated  (Eq.  3). 
Similar  arguments  can  be  made  for  amplitude  data  (A);  however, 
^medium  >5  determined  by  evaluating  measured-  and  instrument- 
response  ratios.  The  instrument  response  was  calculated  by 
making  a  measurement  on  a  homogenous  phantom  with  known 
optical  properties.  A  siloxane  phantom  cast  from  flexible  sili¬ 
cone  (RTV  615/700,  General  Electric)  and  titanium  diode  in  a 
400-ml  mold  was  used  to  calibrate  the  instrument  response  in 
phase  and  amplitude  measurements.  Optical  properties  of  the 
phantom  were  determined  by  the  two-distance  technique  de¬ 
scribed  by  Fantini  et  al.  (33).  Reference  measurements  were 
made  on  the  phantom  at  the  given  source-detector  separation 
during  the  time  of  the  subject  measurement.  Calibrations  were 
made  at  approximately  10-min  intervals  to  account  for  ampli¬ 
tude  drift. 
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1,897  X  10®  2.037  X  10®  0.34 

1.809  X  10®  2.659  X  10®  0.781 
1.570X  10®  3.049  X  10®  0.74 

Given  that  there  are  more  wavelengths  than  principle  chro¬ 
mophores,  we  have  four  equations,  and  three  unknowns,  and  no 
general  solution  for  c,.  The  chromophore  concentrations  are 
determined  by  using  a  least-squares  solution  to  Ec  =  /Xg.  In 
matrix  representation,  the  chromophore  concentration  is  given 
by  using  c  =  EY^  E^^x.^,  where  and  E~^  denote  the 
transpose  and  inverse  of  the  matrix  E,  respectively. 


[Hb] 
[HbO^l 
.  [H2O]  J 


[2] 


Calibration.  A  measurement  of  the  phase  shift  of  intensity- 
modulated  light  at  source  wavelength  A  propagating  through  a 
turbid  medium  at  a  given  source-detector  separation  can  be 
expressed  as  follows: 

^meas  '^medium  T  '^instrument*  [2] 


Results  and  Discussion 

The  results  of  raw  data  and  simultaneous  fits  to  the  P\  approx¬ 
imation  for  Pre  subject  1  and  Post  subject  1  are  shown  in  Fig.  2. 
Clear  differences  are  detected  between  the  two  subjects  in  both 
amplitude  and  phase  data. 

To  establish  the  range  of  optical  property  values  characteristic 
of  normal  menopausal  states,  we  examined  nine  subjects  (six 
premenopausal)  not  receiving  exogenous  hormone  therapy.  Fig. 
3  shows  a  scatter  plot  of  jxa  versus  p,s'  at  all  four  wavelengths  for 
each  subject.  For  comparative  purposes,  data  from  the  right 
upper-outer  quadrant  measurements  are  shown.  Glandular  tis¬ 
sue  is  concentrated  in  the  upper-outer  quadrant  of  the  breast, 
and  consequently,  more  lesions  appear  in  this  region  than  any 
other  (34).  Fig.  3  shows  that,  in  general,  postmenopausal  women 
have  substantially  lower  jxa  and  /Xs'  values  compared  with 
premenopausal  women.  Average  premenopausal  p,a  values  at 
each  wavelength  (0.0048-0.015  mm"^)  are  2.3-  to  3-fold  higher 
than  postmenopausal  absorption  (average  p,a  =  0.0016-0.0064 
mm"’^).  Furthermore,  all  postmenopausal  pa  values  are  less  than 
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Absorption  vs.  Scattering 
for  Post  and  Pre  subjects 
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Fig.  3.  Iia  vs.  for  six  pre-  and  three  postmenopausal  subjects  at  all 
wavelengths.  Values  are  calculated  from  best  diffusion-model  fits  to  phase 
and  amplitude  data 
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Fig.  4.  fis  vs.  wavelength  for  three  subjects  of  varied  hormonal  and  meno¬ 
pausal  status.  Some  values  at  980  nm  represent  values  extrapolated  from 
linear  regression  analysis  of  the  data. 


0.007  mm~’.  Premenopausal  women  have  16-22%  higher  .scat¬ 
tering  values  at  each  wavelength  (average  /Xs'  =  0.83-1.1  mm"') 
than  postmenopausal  subjects  (average  /Xs'  =  0.67-0.86  mm"'). 
Because  measurement  uncertainties  are  less  than  5%  of  the 
optical  property  value  (29),  these  valves  constitute  significant 
and  dramatic  differences  between  menopausal  states. 

A  deeper  understanding  of  the  underlying  physiological  rea¬ 
sons  for  optical  property  contrast  can  be  obtained  by  examining 
jXs'  spectra  (i.e.,  /x/  vs.  A)  for  three  individuals.  Results  shown  in 
Fig.  4  highlight  the  impact  of  estrogen  on  breast  structure.  The 
29-year-old  premenopausal  subject  displays  p,s'  values  that  are 
consistently  greater  than  the  67-year-old  postmenopausal  sub¬ 
ject’s  values.  Interestingly,  /Xs'  values  for  the  52-year-old  post¬ 
menopausal  subject  receiving  estrogen-only  HRT  fall  directly 
between  premenopausal  and  postmenopausal  data.  The  elevated 
HRT  valves  may  be  caused  by  the  fact  that  estrogen  increases  the 
rate  of  mitosis  within  ductal  tissue.  Thus,  optical  property 
differences  between  HRT  and  postmenopausal  subjects  could  be 
a  consequence  of  epithelial  tissue  proliferation.  Scattering  con¬ 
trast  between  premenopausal  and  postmenopausal  subjects  is 
probably  caused  by  a  combined  effect  of  the  loss  of  glandular 
epithelium  as  well  as  extracellular  matrix  remodeling.  The  more 
gradual  pis'  vs.  A  slope  shown  for  the  postmenopausal  subjects 
reflects  large  particle  scattering  consistent  with  a  high  percent¬ 
age  of  fatty  parenchyma  found  in  this  age  group.  In  younger 
women,  the  relatively  steep  /Xs'  vs.  A  slope  is  likely  to  be 
influenced  by  the  presence  of  extracellular  collagen  in  addition 
to  cellular/epithelial  factors. 

These  concepts  are  substantiated  further  by  Fig.  5,  which 
illustrates  an  inverse  correlation  between  jx/  and  age  for  subjects 
over  50.  Whether  the  trend  indicating  increased  scattering  for 
HRT  women  is  caused  by  exogenous  hormone  use  is  not  yet 
clear.  Premenopausal  values  seem  to  be  independent  of  age, 
most  likely  because  of  the  timing  of  the  measurements.  The  six 
premenopausal  subjects  were  measured  at  different  points  of 
their  menstrual  cycles,  which  can  significantly  affect  optical 
properties. 


Average  tissue  hemoglobin  concentrations  (HbO:,  deoxy,  and 
total)  and  tissue  water  concentration  are  shown  in  Fig.  6.  The 
concentration  of  blood  vessels  and  blood  flow,  indicated  by  total 
hemoglobin  concentration  (THC)  (THC  =  [Hb]  -I-  [HbO^]),  is 
greatest  in  premenopausal  breast  because  of  the  high  vascular 
demands  of  the  tissue.  Blood  vessel  density  and  blood  flow 
decrease  as  a  result  of  menopausal  involution  of  glandular  tissue. 
This  effect  can  be  detected  from  the  significant  total  hemoglobin 
concentration  difference  between  premenopausal  and  post¬ 
menopausal  subject  groups. 


^  ^  Scattering  vs.  Age 
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Fig.  5.  at  674  nm  vs.  age  for  all  subjects.  Pre  6  (age  18)  is  represented  as 
an  average  of  two  points  corresponding  to  two  separate  measurement  dates. 
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Fig.  6.  Mean  Hb  concentrations  [^r,M,  Hb02,  deoxy,  and  total  (THC)]  and 
water  (H2O)  concentration  relative  to  pure  water  (%)  for  each  subject  group. 
Values  are  determined  from  wavelength-dependent  absorption  values  at  674, 
803,  849,  and  956  or  980  nm.  Error  bars  represent  the  normalized  standard 
deviation  to  the  mean  for  six  premenopausal,  five  HRT,  and  three  postmeno¬ 
pausal  subjects.  Confidence  values  are  for  Hb,  Hb02,  THC,  and  H2O. 


Women  using  HRT  exhibit  Hb  concentration  values  that  are 
intermediate  between  premenopausal  and  postmenopausal  sub¬ 
jects.  Thus,  FDPM  measurements  may  be  sensitive  to  subtle 
HRT  effects,  including  elevations  in  blood  flow,  fibroglandular 
volume,  and  epithelial  cell  proliferation  (35-38).  The  general 
decrease  in  total  Hb  concentration  with  age  after  50  (Fig.  7) 
suggests  that  it  is  difficult  to  determine  from  limited  existing  data 
whether  the  measured  HRT  Hb  effect  is  real  or  simply  a 
consequence  of  age-related  optical  property  changes. 
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Fig.  7.  Total  Hb  concentration  vs.  age  for  all  subjects.  Pre  6  (age  18)  is 
represented  as  an  average  of  two  points  corresponding  to  two  separate 
measurement  dates. 


Table  1.  Variation  in  physiological  properties  during  the 
menstrual  cycle 

Parameter  Day  14  Day  25 

Hb,  pM  3.07  ±  0.04  4.72  ±  0.05 

Hb02,  pM  7.56  ±  0.04  1 1 .04  ±  0.06 

Total  Hb,pM  10.63  +  0.05  15.76  +  0.08 

5,02,  %  71.12  ±  0.53  70.06  ±0.53 

Water,  %  16.49  ±0.11  21.12  ±0.15 

Calculated  physiological  properties,  Hb  concentration  (pM,  Hb02,  deoxy-, 
and  total  Hb),  blood  oxygen  saturation  (5,02),  and  water  concentration 
relative  to  pure  water  (%)  for  a  premenopausal  subject  at  Days  14  and  25  of 
a  28-day  menstrual  cycle. 

Fig.  6  data  show  that  the  average  water  concentration  for 
premenopausal  women  (26  ±  6%)  is  more  than  a  factor  of  two 
greater  than  the  mean  tissue  water  concentration  for  postmeno¬ 
pausal  women  (11  ±  2%).  These  differences  reflect  the  high 
water  content  of  epithelial  connective-tissue  compartments  in 
premenopausal  tissue,  whereas  the  postmenopausal  breast  is 
dominated  by  low  water-content  adipose.  Water  concentration 
in  HRT  subjects  is  slightly  greater  than  the  postmenopausal 
group’s  concentration  (15  ±  7%).  Differences  may  be  caused  by 
hormone-induced  accumulation  of  fluids,  a  commonly  occurring 
side  effect  from  HRT  (38). 

Investigating  positional  variations  in  optical  properties  be¬ 
tween  pre-  (Pre  1)  and  postmenopausal  (Post  1)  subjects  reveals 
a  higher  degree  of  variation  in  scattering  for  post-  («“8%  for  all 
wavelengths)  vs.  premenopausal  women  («=4%).  These  values 
were  calculated  by  normalizing  the  standard  deviation  of  the 
values  from  measurements  in  each  of  the  four  quadrants  (n  =  4) 
to  the  mean  (data  not  shown).  This  difference  may  represent  the 
nonuniform  glandular  involution  of  breast  tissue  that  accompa¬ 
nies  menopause  and  results  in  palpable  differences  in  density  in 
the  breast  tissue  of  postmenopausal  subjects  (39).  The  positional 
variation  in  the  p-a  is  greater  for  wavelengths  corresponding  to 
hemoglobin  absorption  (674-849  nm)  for  the  postmenopausal 
subjects  (23-33%  vs.  17-20%  for  Pre).  However,  wavelengths 
corresponding  to  fat  and  water  concentrations  (894, 947,  and  956 
nm)  show  a  higher  variability  in  premenopausal  (13-16%)  vs. 
only  7-8%  in  postmenopausal  subjects. 

Figs.  3-7  strongly  suggest  that  physiological  changes  caused  by 
hormonal  fluctuations  that  occur  over  a  period  of  many  years  can 
be  detected  and  quantified.  Previous  studies  (40)  show  these 
effects  also  are  detectable  within  the  menstrual  cycle  of  pre¬ 
menopausal  women.  To  test  FDPM  sensitivity  to  menstrual- 
cycle  variations,  we  examined  a  premenopausal  subject  during 
ovulation  and  before  the  onset  of  menses.  Table  1  data  summa¬ 
rizing  these  results  show  that  /j-a  and  /j-s'  are  higher  for  the  latter 
part  of  the  cycle  (Day  25)  than  in  mid-cycle  (Day  14).  Absorption 
differences  correspond  to  the  calculated  physiological  proper¬ 
ties.  Our  results  show  a  48.3%  increase  in  Hb  and  28.1%  increase 
in  water  concentration  during  the  luteal  phase,  changes  that  are 
consistent  with  the  physiological  effects  caused  by  ovarian 
hormone  fluctuations  during  the  menstrual  cycle.  After  ovula¬ 
tion,  blood  flow  to  the  breast  can  increase  by  up  to  50%;  there 
is  an  increase  in  breast  volume;  and  parenchymal  water  content 
changes  by  an  average  of  25%  during  the  latter  half  cycle  (24, 25). 
Changes  in  scattering  were  3-5%  and  not  statistically  significant. 

Conclusions 

A  hand-held  photon  migration  probe  has  been  developed  that 
can  detect  significant  differences  in  absorption  and  scattering 
properties  occurring  with  changes  in  age  and  menstrual-cycle 
phase.  The  sensitivity  of  quantitative  NIR  spectroscopy  to  breast 
biology  is  unique  among  radiological  methods.  Consequently, 
optical  techniques  may  eventually  provide  a  practical,  noninva- 
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sive  means  for  enhancing  the  accuracy  of  tumor  diagnostics  in 
premenopausal  and  perimenopausal  subjects  and  for  increasing 
our  understanding  breast  pathophysiology. 
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Broadband  absorption  spectroscopy  in  turbid 
media  by  combined  frequency-domain  and 
steady-state  methods 


Frederic  Bevilacqua,  Andrew  J.  Berger,  Albert  E.  Cerussi,  Dorota  Jakubowski,  and 
Bruce  J.  Tromberg 


A  technique  for  measuring  broadband  near-infrared  absorption  spectra  of  turbid  media  that  uses  a 
combination  of  frequency-domain  (FD)  and  steady-state  (SS)  reflectance  methods  is  presented.  Most  of 
the  wavelength  coverage  is  provided  by  a  white-light  SS  measurement,  whereas  the  FD  data  are  acquired 
at  a  few  selected  wavelengths.  Coefficients  of  absorption  (|jl„)  and  reduced  scattering  (p.')  derived  from 
the  FD  data  are  used  to  calibrate  the  intensity  of  the  SS  measurements  and  to  estimate  p'  at  all 
wavelengths  in  the  spectral  window  of  interest.  After  these  steps  are  performed,  one  can  determine  p„ 
by  comparing  the  SS  reflectance  values  with  the  predictions  of  diffusion  theory,  wavelength  by  wave¬ 
length.  Absorption  spectra  of  a  turbid  phantom  and  of  human  breast  tissue  in  vivo,  derived  with  the 
combined  SSFD  technique,  agree  well  with  expected  reference  values.  All  measurements  can  be  per¬ 
formed  at  a  single  source-detector  separation  distance,  reducing  the  variations  in  sampling  volume  that 
exist  in  multidistance  methods.  The  technique  uses  relatively  inexpensive  light  sources  and  detectors 
and  is  easily  implemented  on  an  existing  multiwavelength  FD  system.  ©  2000  Optical  Society  of 
America 

OCIS  codes:  170.1470,  170.4090,  170.5280,  170.7050. 


1.  Introduction 

Reflectance  spectroscopy  is  a  technique  for  character¬ 
izing  turbid  media  that  has  become  widely  used  in 
medical  diagnostics.  In  many  cases  the  quantifica¬ 
tion  of  chromophore  concentrations  is  desired,  cmd 
this  requires  the  ability  to  separate  the  effects  of 
absorption  from  those  of  scattering.  Fundamen¬ 
tally,  the  coefficients  of  absorption  |ji,„  and  of  reduced 
scattering  can  be  determined  by  a  series  of  reflec¬ 
tance  measurements  performed  in  one  of  three  do¬ 
mains,  namely,  time^-^  (with  a  fast  pulse  of  light), 
frequency^-'^  (with  a  sinusoidally  modulated  source  of 
light),  and  steady  state®-^®  (with  a  source  of  constant 
intensity  but  multiple  detectors  at  different  distanc¬ 
es).  Unsurprisingly,  these  three  techniques  have 
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different  merits  and  limitations.  Spatially  re¬ 
solved  steady-state  techniques  are  relatively  inex¬ 
pensive  and  are  more  readily  suited  for  the 
determination  of  ijl„  and  [x'  over  large,  continuous 
ranges  of  wavelengths  than  are  the  other  methods. 
However,  the  steady-state  approach  works  best 
when  measurements  are  performed  with  a  combi¬ 
nation  of  short  (~1  transport  mean  free  path)  and 
long  (many  transport  mean  free  paths)  source- 
detector  separations. Ideally,  the  optical  proper¬ 
ties  of  the  sample  should  not  vary  over  the  ranges  of 
volumes  probed  by  the  various  measurements. 
The  larger  the  spread  of  distances  probed,  the  more 
likely  that  heterogeneities,  such  as  those  found  in 
biological  tissue,  will  distort  the  data  from  the  pre¬ 
dictions  of  the  model.  One  approach  to  limiting 
this  effect,  given  that  the  shortest  separations  pro¬ 
vide  great  stability  for  the  calculation  of  (x^,  is  to  use 
relatively  short  (<10-mm)  source- detector  separa¬ 
tions. Inasmuch  as  the  mean  probing  depth 
scales  with  the  source— detector  separation,  with 
this  approach  such  measurements  are  sensitive  to 
superficial  components  (to  depths  of  less  than  5  mm 
for  typical  biological  tissues). 

Time-  and  frequency-domain  techniques  are  well 
suited  for  deeper  (>1  cm  for  biological  tissue)  inves- 
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tigations.  Moreover,  they  can  be  performed  with 
only  one  or  a  few  source- detector  separations,  which 
makes  them  more  robust  for  use  in  studying  hetero¬ 
geneous  samples.  Because  such  techniques  require 
sources  that  can  be  pulsed  or  modulated  rapidly,  cov¬ 
ering  a  large  wavelength  range  requires  a  tunable 
laser  or  an  extensive  collection  of  laser  diodes,  both  of 
which  can  be  expensive,  difficult  to  maintain,  and 
slow  to  cover  the  entire  spectrum.  This  is  an  impor¬ 
tant  drawback,  because,  as  discussed  by  Hull  et 
the  quantification  of  chromophore  concentrations  can 
be  significantly  affected  by  use  of  a  limited  number  of 
wavelengths. 

In  this  paper  we  suggest  a  way  to  use  steady-state 
(SS)  and  frequency-domain  (FD)  reflectance  mea¬ 
surements  in  tandem  to  obtain  broad  wavelength  cov¬ 
erage  with  increased  penetration  depth.  This 
method  is  especially  promising  for  near-infrared 
spectroscopy  of  tissue,  e.g.,  to  characterize  breast 
physiology.  For  such  applications,  the  method  pro¬ 
posed  here  permits  rapid  data  acquisition,  deep  tis¬ 
sue  probing,  and  robust  resolution  of  the 
contributions  from  the  four  major  near-infrared  tis¬ 
sue  absorbers:  oxyhemoglobin,  deoxyhemoglobin, 
water,  and  fat.  The  central  innovations  are  using 
FD-derived  and  values  to  convert  the  SS  mea¬ 
surements  into  units  of  absolute  reflectance  and  us¬ 
ing  the  power-law  wavelength  dependence  of  p'  to 
obtain  interpolated  and  extrapolated  values  at  non¬ 
laser  wavelengths.  FD  measurements  are  made  at 
a  handful  of  diode  laser  wavelengths  spanning  the 
range  of  interest  (650-1000  nm),  whereas  the  SS 
measurements  are  made  continuously  across  the  en¬ 
tire  range.  Unlike  for  spatially  resolved  SS,  how¬ 
ever,  here  only  a  single,  large  source- detector 
separation  is  used,  preferably  the  same  one  as  for 
the  FD  measurements.  The  instrumentation  is 
straightforward  and  particularly  easy  to  add  to  an 
existing  FD  system. 

After  describing  our  method  of  combining  the  SS 
and  FD  methods  (SSFD),  we  test  the  SSFD  technique 
by  using  it  to  measure  the  absorption  spectra  of  tur¬ 
bid  samples.  First,  we  analyze  a  liquid  tissue  phan¬ 
tom  whose  absorption  spectrum  is  loaown  (by  direct 
spectrophotometry  of  the  absorbing  component  be¬ 
fore  mixing)  and  with  which  the  SSFD  result  can  be 
compared.  We  also  measure  locations  on  the  breasts 
of  two  human  female  volunteers,  demonstrating  that 
data  can  be  gathered  in  vivo  and  analyzed  to  provide 
pertinent  physiological  parameters.  Differences 
among  the  measured  breast  spectra  can  be  inter¬ 
preted  in  terms  of  different  relative  levels  of  water, 
fat,  and  hemoglobin  present  in  the  explored  tissue 
volumes.  In  addition,  estimates  of  absolute  concen¬ 
trations  are  comparable  with  those  reported  in  other 
recent  broadband  in  vivo  studies.  Finally,  we  dis¬ 
cuss  possible  reasons  for  imperfect  spectral  fits  and 
compare  concentration  predictions  derived  from 
SSFD  analysis  with  those  derived  from  FD  data 
alone. 
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Fig.  1.  Configuration  of  light  sources,  optical  fibers,  and  detec¬ 
tors:  APD,  avalanche  photodiode.  The  dashed  rectangles  denote 
components  that  belong  to  the  FD  and  SS  systems. 


2.  Experimental  Methods 

A.  Optical  Measurements 

Figure  1  shows  the  experimental  arrangement  for 
SSFD  measurements.  In  all  cases,  light  is  delivered 
via  optical  fiber  to  the  surface  of  the  sample  and 
collected  at  some  distance  p  away.  For  the  liquid 
phantom  measurements,  p  was  15.5  mm;  for  the 
breast,  21.5  (FD)  and  24  (SS)  mm  (the  slight  differ¬ 
ence  was  due  to  instrumental  limitations;  a  future 
instrument  will  utilize  identical  distances).  In  FD 
mode  (upper  dashed  rectangle),  the  light  arrives  se¬ 
quentially  from  one  of  seven  amplitude-modulated 
diode  lasers  (672,  800,  806,  852,  896,  913,  and  978 
nm,  all  with  output  powers  of  <20  mW  at  the  sample) 
and  is  detected  by  an  avalanche  photodiode  unit 
(Hamamatsu  C556P-56045-03)  that  amplifies  the  ac 
component  of  the  signal.  A  network  analyzer 
(Hewlett-Packard  8753C)  delivers  251  modulation 
frequencies  from  100  to  700  MHz  and  measures 
phase  and  modulation  amplitudes  of  the  photon  in¬ 
tensity  signal,  as  described  elsewhere.'^-^®  In  SS 
mode  (lower  dashed  rectangle),  light  comes  from  a 
150-W  halogen  lamp  (Fiber-Lite)  and  is  analyzed  by  a 
fiber-coupled  spectrograph  (Ocean  Optics  S2000) 
with  a  linear  CCD  detector  from  525  to  1155  nm,  with 
the  useful  range  for  our  experiments  being  650-1000 
nm.  The  spectrograph  records  a  total  of  2048  points 
(0.35  nm/pixel),  and  the  spectral  resolution  is  5  nm 
(full  width  at  half-maximum).  Light  is  delivered  to 
the  sample  through  a  bundle  of  four  fibers  (bundle 
diameter,  600  p.m)  and  collected  with  a  single  fiber  of 
1-mm  diameter.  We  measure  the  spectrum  of  the 
light  source  separately  by  inserting  the  source  and 
detector  fibers  into  different  ports  of  an  integrating 
sphere  (Labsphere,  IS-040-SF).  Relative  reflectance 
is  calculated  to  be  the  sample  spectrum  divided  by  the 
source  spectrum  (note  that  both  measurements  use 
the  same  delivery  fiber,  collection  fiber,  and  detector 
apparatus).  Total  acquisition  time  per  sample  for 
SSFD  measurements  is  of  the  order  of  40  s  (30  s  for 
FD  and  10  s  for  SS).  We  calculated  p.a(X)  according 
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to  the  methods  of  Section  3  with  an  in-house  Matlab 
(The  MathWorks,  Inc.)  code,  making  use  of  the  opti¬ 
mization  toolbox. 

B.  Samples 

The  phantom  contained  17  mg  of  green  dye  (naph- 
thol)  and  50  mL  of  an  aqueous  scattering  suspension 
(Intralipid-20%,  Pharmacia,  Inc.)  dissolved  in  950 
mL  of  water.  Before  the  addition  of  Intralipid,  the 
absorbance  spectrum  of  the  dye-water  solution  was 
measured  in  a  1-cm  plastic  cuvette  by  a  spectropho¬ 
tometer  (Beckman  Instruments,  DU630)  with  water 
as  a  reference.  Optical  measurements  were  per¬ 
formed  as  indicated  in  Fig.  1,  with  source  and  detec¬ 
tor  fibers  placed  at  the  surfaces  of  the  samples. 

In  vivo  measurements  were  performed  on  the  left 
breasts  of  two  supine  female  volunteers,  aged  37  and 
21  years.  Data  were  gathered  from  two  regions  on 
one  volunteer  (a  region  on  the  areolar  border  and  a 
region  of  inner  breast,  i.e.,  close  to  the  middle  of  the 
chest)  and  from  one  region  (inner  breast  only)  on  the 
other  volimteer,  jfielding  a  total  of  three  samples.  In 
this  modality,  the  source  light  was  again  delivered  by 
optical  fiber  but  the  FD  detector  was  placed  directly 
against  the  tissue,  without  a  collection  fiber.  Fiber 
and  detector  were  bimdled  into  a  single  hand-held 
device  that  was  placed  gently  against  the  breast.  SS 
reflectance  was  measured  subsequently,  at  the  exact 
same  location,  in  the  two-fiber  mode  described  above. 
All  procedures  were  approved  by  the  Institutional 
Review  Board  of  the  University  of  California,  Irvine 
(study  95-563). 


where  to  is  the  modulation  frequency  in  radians  per 
second.  The  steady-state  solution  is  simply  the  limit 
when  0)  =  0.  In  this  case  the  solution  for  i?  is  a  real 
number  {^^eai  =  [3p-„(p-a  +  and  =  0}. 

When  (a  >  0,  R  becomes  a  complex  number  A 
exp(- jc()),  where  A  is  the  modulation  amplitude  and  c|) 
is  the  phase  shift  relative  to  the  source.  As  we  ex¬ 
plain  in  Subsection  3.B,  these  quantities  can  be  ob¬ 
tained  from  the  FD  measurement. 

In  treating  reflectance  problems,  we  model  the 
sample-air  interface,  using  an  extrapolated  bound¬ 
ary  conditional  in  which  the  fluence  is  set  to  zero  at  a 
distance  Zj,  =  2D(1  -I-  i?eg)/(l  -  Rgg)  above  the  sam¬ 
ple  (Rgjj-  depends  on  the  refractive  mismatch  and 
equals  0.493  for  tissue  of  n  =  1.4  and  air  of  n  =  1.0). 
The  method  of  images  is  employed,  with  an  isotropic 
point  source  at  a  depth  =  (p,„  -I-  |jLs)“^  contributing 
a  signal  and  a  negative  image  point  at  a  height 
Zg  -t  Z/,  above  the  extrapolated  boundary  contributing 
Sj.  The  result  is  the  fluence  d>  at  any  point  in  the 
sample.  The  detected  signal  along  the  boundary  is 
then  written  as  a  combination  of  terms  proportional 
to  the  fluence  and  to  its  flux  normal  to  the  surface: 

R  =  cid)  -  CjDVO)  •  i-z),  (1) 

where  values  for  the  constants  and  C2  are  deter¬ 
mined  by  the  refractive-index  mismatch  between  the 
two  media^i  (for  tissue  of  n  =  1.40  and  air  of  n  =  1, 
for  these  constants  the  values  0.118  and  0.306,  re¬ 
spectively,  are  assumed),  -2  is  a  unit  vector  pointing 
normally  upward  out  of  the  sample,  and 


3.  Background  Theory 

A.  Diffusion  Model 

The  measured  reflectance  signal  R  can  be  predicted 
theoretically  by  use  of  the  diffusion  approximation  to 
the  radiative  transfer  equation,  as  many  groups  of 
researchers  have  discussed.'^-i'^-^i  In  this  theoretical 
framework,  the  reflectance  is  a  function  of  the  opticEd 
properties  of  the  medium,  defined  by  absorption  co¬ 
efficient  |x„,  reduced  scattering  coefficient  p,',  and  in¬ 
dex  of  refraction  n.  The  diffusion  approximation  is 
valid  for  large  source- detector  separation  r  [r  > 
10(p„  -I-  Ps)~^]  and  high  (reduced)  albedo  [p.y(p.a  + 
Ps)  >  0.95].  The  SS  and  FD  cases  can  be  described 
by  a  single  formalism  in  which  the  solution  for  the 
reflectance  is  built  from  the  Green’s  function  for  the 
diffusion  equation,  i.e.,  the  fluence  [W/cm^]  that  is 
due  to  an  isotropic  point  source  in  an  infinite,  homo¬ 
geneous  medium.  This  function  takes  the  form 
exp{-kr)/{Dr),  where  D  =  [3(pa  -I-  Pg)]”^,  the  com¬ 
plex  wave  number  is  ^  =  ^^eai  + 
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where  P  is  the  incident  power  and  (rj)  is  the  dis¬ 
tance  from  the  source  (image)  to  the  detector. 


B.  Frequency-Domain  Fitting 

Like  the  theory  of  Subsection  3. A,  our  FD  fitting  pro¬ 
cess,  which  provides  p„  and  p'  values  at  a  few  wave¬ 
lengths,  has  been  described  at  length  elsewhere,'^-!® 
and  a  brief  description  is  provided  here  as  back¬ 
ground. 

Each  FD  measurement  contains  instrumental  ar¬ 
tifacts;  i.e.,  the  measured  reflectance  is  actually 


R  =  CoA  exp[-i(c|)  -h  4)o)],  (3) 
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where  and  <|)o  £u-e  sample-independent  instrumen¬ 
tal  constants.  A  and  cj)  are,  respectively,  the  modu¬ 
lation  amplitude  and  phase  of  the  FD  reflectance,  as 
defined  in  Subsection  3. A. 

For  calibration,  we  gather  FD  data  from  a  prepared 
sample  whose  p„  and  p^  values  are  known  a  priori 
from  a  set  of  two-distance  FD  measurements.'^ 
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Frequency-dependent  values  of  (j)^  and  C„  are  calcu¬ 
lated  from  the  discrepancies  between  measured  [Eq. 
(3)]  and  predicted  [Eq.  (1)]  phase  and  modulation 
amplitude,  thus  calibrating  our  future  measure¬ 
ments. 

With  the  instrumental  constants  thus  determined, 
the  502  data  points  per  sample  (251  for  both  phase 
and  amplitude  as  a  function  of  the  FD  modulation 
frequency  w)  depend  on  two  unknowns,  p.„  and  p,'. 
We  select  the  best  and  p'.  values  to  fit  the  predic¬ 
tions  of  Eq.  (1)  to  the  data,  using  the  iterative,  non¬ 
linear,  least-squares  method  of  Levenberg  and 
Marquardt,^^  simultaneously  fitting  the  phase  and 
the  amplitude  imtil  overall  convergence  is  achieved.'' 
As  was  noted  by  various  authors,®’^'  the  reflectance  in 
Eq.  (1)  is  nearly  proportional  to  either  the  flux  or  fhe 
fluence  term  alone  (i.e.,  the  ratio  of  the  two  is  inde¬ 
pendent  of  p„  and  p')  for  the  large  source- detector 
separations  that  we  employ  [>10(p<j  +  p^)"^].  We 
therefore  use  only  the  fluence  term  <I>  from  Eq.  (1), 
absorbing  the  additional  proportionality  factor  into 
the  instrumental  constant  C^. 

4.  Determination  of  Broadband  p^  Spectrum 

Section  3  provided  background  information  regard¬ 
ing  diffusion  theory  models  and  a  method  for  extract¬ 
ing  p„  and  Ps  values  from  FD  measurements  at 
specific  wavelengths.  We  now  describe  how  FD  and 
SS  measurements  can  be  combined  to  5deld  quanti¬ 
tative,  broadband  p^  spectra. 

Our  goal  is  to  compute  p^  at  each  wavelength, 
given  the  measured  SS  reflectance.  However, 
single-distance  SS  reflectance  cannot  itself  provide 
p^:  R  depends  on  p^  as  well  as  on  p^,  so  one  mea¬ 
surement  cannot  provide  a  unique  determination  of 
either  parameter.  In  addition,  the  instrumental 
constant  is  not  known.  The  single  SS  measure¬ 
ment  at  each  wavelength  therefore  needs  to  be  sup¬ 
plemented  with  both  a  means  of  establishing  the 
absolute  reflectance  intensity  and  an  additional  piece 
of  information  about  p^  or  p^  (or  a  combination  of  the 
two). 

A.  Use  of  Frequency-Domain  Data  to  Provide  Additional 
Information 

The  FD  system  can  provide  both  pieces  of  necessary 
information  at  each  SS  wavelength.  At  first  glance 
this  is  surprising,  as  the  FD  system  operates  at  only 
seven  wavelengths,  whereas  the  SS  system  covers 
450  nm  with  continuous  5-nm  resolution.  It  would 
seem  then  that  the  needed  values  of  p'  and  C„  could 
be  determined  only  at  these  seven  wavelengths,  and 
strictly  speaking  such  is  indeed  the  case.  Signifi¬ 
cantly,  however,  the  wavelength  dependence  of  both 
Ps  and  Cq  is  smooth  and  predictable  in  shape,  as  is 
discussed  below.  It  is  this  fact  that  enables  us  to  use 
the  discrete  FD  information  to  supplement  our  broad¬ 
band  SS  measurements  at  all  desired  wavelengths. 

B.  Amplitude  Calibration  of  the  Reflectance  Spectrum 

The  FD  measurements  permit  calculation  of  the  in¬ 
strumental  factor  Cg  at  all  wavelengths  because  we 


expect  no  wavelength  dependence  at  all:  The  sam¬ 
ple  and  source  spectra  are  measured  with  the  same 
delivery  and  collection  system  and  only  a  few  seconds 
apart,  so  there  should  be  negligible  wavelength- 
dependent  artifacts  in  the  ratio  of  the  two.  There¬ 
fore  the  task  reduces  to  that  of  calculating  at  a 
single  wavelength.  This  is  readily  done  at  any  of  the 
FD  wavelengths,  as  p^  and  p^  are  both  known  and 
absolute  reflectance  is  a  function  only  of  these  vari¬ 
ables  [cf  Eqs.  (1)  and  (2)  in  the  w  =  0  limit].  Using 
all  the  FD  data  to  increase  robustness,  we  can  calcu¬ 
late  the  value  of  that  scales  the  measured  SS 
reflectances  to  match  the  predicted  reflectances  as 
closely  as  possible  (in  the  least-squares  sense). 
Once  this  is  done,  the  scale  factor  is  known  for  the 
entire  spectrum.  Examples  of  this  scaling  are  given 
below. 

C.  Determination  of  P3(\) 

The  FD  values  of  p'  also  allow  us  to  obtain  informa¬ 
tion  about  Ps  across  the  entire  wavelength  range. 
As  numerous  groups  of  researchers  have 
observed,23-25  the  particle  size  distribution  of  scatter- 
ers  (0.1-10  pm)  in  many  biological  media  and  phan¬ 
toms  tends  to  have  smooth  wavelength  dependence 
over  the  range  650-1000  nm,  which  is  well  described 
by  a  power  function  of  the  form 

p((X)=A\'^  (4) 

(the  parameter  A  here  has  no  relation  to  the  modu¬ 
lation  amplitude).  As  a  result,  measuring  several  p^ 
values  allows  us  to  fit  those  values  to  a  simple  func¬ 
tion  of  wavelength  and  to  obtain  good  estimates  of  p' 
at  all  other  wavelengths  needed. 

D.  Iterative  Solving  for  the  Paf\)  Spectrum 

With  p'(X)  calculated  and  the  reflectance  correctly 
scaled,  the  equation  for  steady-state  reflectance  ([Eq. 
(1)]  contains  only  one  unknown,  p„.  Because  it  is 
difficult  to  obtain  an  analytical  solution  for  p^  from 
this  equation,  numerical  solving  methods  were  em¬ 
ployed.  Proceeding  one  wavelength  at  a  time,  we 
used  the  Matlab  nonlinear  zero-finding  function  fzero 
to  choose  Pq  such  that  Rmeas  ~  •^thy(|J‘-a)  =  where 
/?meas  the  Calibrated  SS  reflectance  measurement 
and  /?thy(P'a)  is  the  theoretical  reflectance  predicted 
by  Eqs.  (1)  and  (2)  for  a  given  trial  value  of  p^.  Solv¬ 
ing  for  the  entire  p^  spectrum  required  —10  s  on  a 
personal  computer. 

5.  Results 

Figure  2  displays  the  SS  reflectance  measured  from 
the  phantom  and  the  predicted  absolute  reflectance 
calculated  from  diffusion  theory  based  on  the  FD 
measurements  of  p^  and  p^.  Because  of  strong  ab¬ 
sorption  at  wavelengths  above  950  nm,  the  FD  mea¬ 
surement  at  978  nm  was  noisy,  which  in  turn  made 
the  calculation  of  SS  data  at  wavelengths  above  950 
unstable.  Data  are  therefore  presented  for  650-950 
nm  for  the  phantom  only.  As  anticipated,  all  the 
measurements  differ  from  the  predictions  by  essen- 
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Fig.  2.  SS  reflectance  spectrum  acquired  from  a  dye-Intralipid 
phantom,  scaled  to  fit  the  discrete  reflectance  values  (filled  circles) 
predicted  by  Eqs.  (1)  smd  (2)  with  FD-derived  values  of  p.'  and 
This  scaling  causes  the  entire  spectrum  to  be  converted  into  abso¬ 
lute  reflectance  units;  see  text  for  a  discussion.  The  error  in  FD 
reflectance  is  estimated  to  be  ±3  x  10^’/mm~^;  in  SS  reflectance, 
±1  X  10“Vnim“^. 


tially  the  same  scale  factor.  The  entire  SS  spectrum 
is  thus  converted  into  absolute  reflectance  units. 
Similar  results  were  observed  for  the  three  breast 
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Fig.  4.  Determination  of  the  broadbcmd  p(  spectrum  for  the  dye- 
Intralipid  phantom.  Open  circles,  discrete  values  measured  by 
the  FD  technique;  solid  curve,  best  power-law  fit  to  Eq.  (4).  Fit¬ 
ting  parameters  are  A  =  2200  and  B  =  -0.82.  Error  bars  on  p( 
are  3%.’ 

of  2  over  the  entire  range,  so  the  wavelength  depen¬ 
dence  is  weak  in  addition  to  being  predictable. 

Figure  6  summarizes  the  absorption  values  calcu¬ 
lated  for  the  dye-Intralipid  phantom.  The  thicker 
curve  shows  the  values  that  we  fotmd  make  Eq.  (1) 
reproduce  the  measured  reflectance  data  as  closely  as 


samples,  as  shown  in  Fig.  3.  From  top  to  bottom, 
sample  1  is  the  inner  breast  region  of  the  37-year-old 
volunteer  and  samples  2  and  3  are  from  the  areolar 
border  and  the  inner  breast,  respectively,  of  the  21- 
year-old  volunteer. 

Figure  4  shows  the  least-squares  power-law  flt  to 
the  iXj  values  measured  by  the  FD  on  the  phantom. 
As  expected,  the  wavelength  dependence  is  fairly 
smooth  and  is  easily  described  by  the  flt.  Corre¬ 
sponding  fits  for  the  three  breast  samples  appear  in 
Fig.  5.  Note  that  the  values  change  by  only  a  factor 


possible,  using  the  fzero  algorithm  as  described 
above.  The  thinner  curve  represents  a  linear  least- 
squares  fit  to  the  thicker  curve  based  on  reference 
spectra  of  pure  water  (from  Kou  et  al.^^)  and  pure  dye 
(measured  by  a  spectrophotometer  as  described 
above).  As  the  figure  shows,  the  two-component  fit 
accurately  models  the  measured  spectrum  across  the 
entire  spectral  range.  The  least-squares  fitting  co¬ 
efficient  for  the  dye  spectrum  is  1.00.  Also  shown 
are  the  values  of  derived  by  FD  analysis  alone, 
demonstrating  an  agreement  between  the  FD  and 
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Fig.  3.  SS  reflectance  spectra  acquired  from  three  locations  in 
breast  tissue  of  female  volimteers.  The  spectra  have  been  scaled 
to  fit  reflectance  values  (filled  circles)  calculated  from  FD  data  in 
the  same  manner  as  for  Fig.  2.  Errors  are  the  same  as  for  Fig.  2. 
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Fig.  5.  Determination  of  broadband  p,'  spectra  for  the  three 
breast  samples.  Circles,  the  discrete  p'  values  measured  by  the 
FD  technique;  solid  curves,  best  power-law  fits  to  Eq.  (4).  Fitting 
parameters  (A,  B)  from  top  to  bottom  are  (240,  -0.86),  (2700,  - 12), 
and  (250,  —0.85).  Error  bars  on  p'  are  3%.'^ 
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Fig.  6.  Comparison  of  |x„  values  generated  by  FD  (open  circles) 
and  by  SSFD  (thicker  curve)  for  the  dye-Intralipid  phantom. 
Also  shown  is  the  best  fit  (thinner  curve)  to  the  SSFD  data  by  use 
of  the  spectra  of  naphthol  (measured  with  a  spectrophotometer) 
and  of  water  (from  Kou  et  Open  circles,  discrete  p.^  values 

from  FD  alone.  Error  bars,  ±0.0005  mm“'  or  5%,  whichever  is 
larger. 


SSFD  methods  at  these  wavelengths,  as  would  be 
expected  because  the  FD  data  have  been  used  to  cal¬ 
ibrate  the  SS  data. 

Corresponding  plots  that  show  absorption  spectra 
of  the  three  breast  samples  appear  in  Figs.  7-9.  As 
before,  the  thicker  curve  is  the  experimental  data,  the 
thinner  curve  is  a  full-spectrum  fit,  and  the  circles  are 
FD  values.  In  these  cases,  however,  the  fit  is  built 
from  published  spectra  (see  Fig.  10)  of  oxyhemoglo¬ 
bin,  deoxyhemoglobin,^'^  water, ^6  and  fatj^s  which  are 
commonly  regarded  as  the  four  major  absorbers  in 
breast  tissue  in  the  650-1000-nm  range.  Coeffi¬ 
cients  from  the  fit  thus  provide  estimates  of  these 


wavelength  (nm) 

Fig.  7.  p„  predictions  for  the  first  breast  measurement  (inner 
breast;  37-year-old  volunteer).  Thicker  curve,  the  SSFD  data; 
thinner  curve,  least-squares  fit  for  a  superposition  of  Hb,  Hb02, 
water,  and  fat  spectra;  filled  circles,  the  FD  values  (error  bars, 
±0.0005  mm"'  or  5%,  whichever  is  larger);  dotted  curve  the  least- 
squares  fit  when  only  the  FD  values  are  weighted.  Physiological 
parameters  from  the  SSFD  spectral  fit:  total  hemoglobin  concen¬ 
tration,  22  pM;  oxygen  saturation,  73%;  water,  15  g/cm^;  fat,  0.75 
g/cm^.  See  Table  1  for  more  details.  Note  that  the  fat  peak 
between  900  and  950  nm  is  significantly  underfitted  by  the  FD 
calculation. 


Fig.  8.  p„  predictions  for  the  second  breast  measurement  (areolar 
border;  21-year-old  volunteer).  Curves  and  circles  have  the  same 
assignments  as  for  Fig.  7.  Physiological  parameters:  total  he¬ 
moglobin  concentration,  30  pM;  oxygen  saturation,  70%;  water, 
0.51  g/cm’;  fat,  0.42  g/cm®.  See  Table  1  for  more  details. 


components’  concentrations.  The  estimates  are  use¬ 
ful  for  comparing  one  sample  with  another;  no  inde¬ 
pendent  attempt  to  validate  concentrations  has  been 
made  (however,  the  concentrations  are  approxi¬ 
mately  correct  in  magnitude  and,  as  we  note  in  Sec¬ 
tion  6  below,  the  values  are  consistent  with  others 
from  the  recent  literature® For  comparison  pur¬ 
poses,  the  dotted  curves  in  Figs.  7-9  show  the  spec¬ 
tral  fit  obtained  when  only  the  FD  wavelengths  are 
weighted,  as  in  a  customary  FD-only  experiment. 
The  two  spectral  reconstructions  are  clearly  different, 
and  these  different  reconstructions  lead  to  different 
concentration  estimates.  All  the  concentration  pre¬ 
dictions,  both  SSFD  and  FD,  are  listed  in  Table  1, 
along  with  the  percent  deviation  of  the  FD  value  from 
the  SSFD  value.  The  hemoglobin  results  are  re¬ 
ported  in  terms  of  total  hemoglobin  (oxyhemoglobin 
plus  deoxyhemoglobin)  and  oxygen  satimation  (oxy¬ 
hemoglobin  divided  by  the  total). 


wavelength  (nm) 

Fig.  9.  p„  predictions  for  the  third  breast  measurement  (inner 
breast;  21-year-old  volunteer).  Curves  and  circles  have  the  same 
assignments  as  for  Fig.  7.  Physiological  parameters:  total  he¬ 
moglobin  concentration,  19  pM;  oxygen  saturation,  72%;  water,  28 
g/cm®;  fat,  0.56  g/cm®.  See  Table  1  for  more  details. 
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Fig.  10.  Absorption  spectra  of  major  absorbers  in  breast  tissue. 
Upper  left,  oxyhemoglobin,  1  pM;  upper  right,  deoxyhemoglobin,  1 
pM;  lower  left,  water,  1  g/cm^;  lower  right,  fat  (soybean  oil),  0.9 
g/cTC?.  The  hemoglobin  spectra  are  from  Wray  et  the  water 

is  from  Kou  et  al.,^  and  the  fat  is  from  the  doctoral  dissertation  of 


Eker.28 


6.  Discussion 

To  extract  from  a  single  SS  reflectance  spectrum, 
one  has  to  convert  the  data  into  absolute  units  (i.e., 
the  instrumental  factor  must  be  removed)  and 
determine  p-'  values  at  all  wavelengths.  We  have 
demonstrated  that  FD  measurements  at  a  few  wave¬ 
lengths  can  accomplish  these  goals  across  the  entire 
650-1000-nm  spectrum.  This  somewhat  surprising 
result  arises  because  Ps(X.)  follows  a  power  law, 
whereas  Cq  has  no  wavelength  dependence  at  all. 
Two  innovations  that  combine  FD  and  SS  methods 
are  thus  available.  First,  when  the  SS  reflectance  is 
scaled  to  match  FD  predictions  of  absolute  reflec¬ 
tance,  the  entire  reflectance  spectrum  is  automati¬ 
cally  calibrated.  Second,  fitting  the  FD  p'  values  to 
a  wavelength-dependent  function  yields  the  extra  in¬ 
formation  that  one  needs  to  extract  p„  from  diffusion 
theory  [Eqs.  (1)  and  (2)].  We  note  that  the  FD- 
derived  p^  and  p^  values  could  just  as  readily  be 


supplied  by  time-domain  measurements  at  a  single 
source- detector  separation. 

The  interpolation  of  the  reduced  scattering  spectra 
by  use  of  a  fit  is  an  important  component  of  this 
technique.  The  lower  two  parts  of  Fig.  5  show  ex¬ 
cellent  agreement  with  the  power  law.  The  top  part 
of  the  figure  shows  reasonably  good  agreement,  but 
small  discrepancies  (of  the  order  of  5-8%)  can  never¬ 
theless  be  observed.  This  behavior  is  likely  not  to  be 
an  actual  divergence  from  the  power  law  but  is  prob¬ 
ably  due  to  a  coupling  effect  between  the  scattering 
and  absorption  quantification  in  the  FD  fit.  Such 
coupling  appears  when  the  diffusion  model’s  assump¬ 
tions  are  not  sufficiently  fulfilled.  The  tissue  heter¬ 
ogeneity,  i.e.,  the  layered  superficial  structure  and 
the  mix  of  deeper  fat  globules  and  fibrous  tissue, 
could  be  a  possible  source  of  deviation  from  the 
model.  Interestingly,  the  highest  fat  concentration 
was  measured  for  this  breast  location  (see  Table  1). 
As  described  recently  by  Doornbos  et  using  a 
power-law  fit  to  calculate  the  p^  spectrum  can  in  fact 
provide  an  advantage  in  computing  p„.  Indeed,  by 
smoothing  the  spectrum,  the  power-law  fit  reduces 
the  spurious  coupling  between  scattering  and  absorp¬ 
tion  caused  by  inaccurate  modeling. 

Compared  with  discrete-wavelength  FD  measure¬ 
ments,  access  to  a  continuous  absorption  spectnun 
offers  two  important  advantages:  better  chro- 
mophore  identification  and  improved  concentration 
quantification.  The  ability  of  the  method  to  recover 
the  continuous  shape  of  the  true  absorption  spectrum 
is  evidenced  by  the  phantom  measurement.  Figure 
6  shows  that  the  experimental  spectrum  is  accurately 
fitted  by  use  of  known  spectra  of  dye  and  water.  The 
fitting  coefficient  of  1.00  for  the  dye  additionally  dem¬ 
onstrates  the  accurate  recovery  of  chromophore  con¬ 
centrations. 

The  breast  spectra  (Figs.  7-9)  further  illustrate 
the  usefulness  of  the  SSFD  technique  for  clinical 
investigations.  As  expected,  the  fit  that  uses  oxy¬ 
hemoglobin,  deoxyhemoglobin,  fat,  and  water  ac¬ 
counts  for  most  of  the  absorption  in  the  650- 


Table  1.  Results  of  Physiological  Predictions  for  Breast  Samples  1-3° 


Breast  Sample 
Number 

Component 

SSFD  Method 

FD  Method 

Only 

Percent  Difference 

1 

Total  Hb  (|xM) 

22 

25 

13 

O2  saturation  (%) 

73 

78 

7 

Water  (g/cm®) 

0.15 

0.14 

-7 

Fat  (g/cm^) 

0.75 

0.44 

-41 

2 

Total  Hb  (p-M) 

30 

30 

0 

O2  saturation  (%) 

70 

79 

13 

Water  (g/cm^) 

0.51 

0.40 

-22 

Fat  (g/cm^) 

0.42 

0.39 

-7 

3 

Total  Hb  (pM) 

19 

19 

0 

O2  saturation  (%) 

72 

81 

13 

Water  (g/cm^) 

0.28 

0.27 

-4 

Fat  (g/cm®) 

0.56 

0.54 

-4 

“Column  5  gives  the  percent  difference  between  the  FD  and  the  SSFD  values,  defined  as  100  x  (FD  -  SSFD)/SSFD. 
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1000-nm  range  in  breast  tissue.  The  quality  of  the 
fit  is  especially  striking  for  wavelengths  longer  than 
800  nm.  In  particular,  the  fat  spectrum  used  in 
the  fit  seems  accurate,  matching  the  main  peak  at 
928  nm  and  the  shoulder  at  895  nm.  These  results 
are  surprisingly  good  considering  that  the  pure-fat 
spectrum  was  measured  in  soybean  oil.^s  These 
measurements  stress  the  important  contribution  of 
fat  to  tissue  absorption,  as  this  chromophore  some¬ 
times  was  neglected  in  previous  research. The 
water  peak  at  976  nm  is  also  well  reproduced  in  the 
measurements.  We  observe  that  the  highest  water 
content  was  measured  for  the  areolar  border  sample 
(sample  2).  Interestingly,  we  note  that  for  both 
dye  and  breast  tissue,  data  fits  to  the  water  spec¬ 
trum  of  Kou  et  al.  were  superior  to  those  of  Hale  and 
Querry,29  particularly  in  the  920-960-nm  regime 
where  the  water  absorption  spectrum  increases 
sharply. 

In  the  650-800-nm  region,  the  major  features  of 
tissue  absorption  are  clearly  due  to  oxyhemoglobin 
and  deoxyhemoglobin.  The  deoxyhemoglobin  peak 
at  760  nm  is  distinctly  visible  in  the  experimental 
data.  Nevertheless,  we  observe  small  but  consistent 
spectral  differences,  of  the  order  of  0.001  mm""^,  be¬ 
tween  fit  and  measurement  in  this  regime  (see  Figs. 
7-9).  These  discrepancies,  revealed  in  the  full  SSFD 
spectrum,  are  not  evident  when  the  FD  data  alone 
are  fitted.  Two  reasons  for  the  imperfect  fits  can  be 
suggested.  First,  the  oxyhemoglobin  and  deoxyhe¬ 
moglobin  absorption  spectra  that  we  used  could  be 
slightly  incorrect.  Indeed,  small  variations  are 
found  among  various  published  spectra,  and  changes 
could  be  also  expected  between  in  vitro  and  in  vivo 
values.  Second,  other,  background,  chromophores 
should  probably  be  taken  into  account.  For  exam¬ 
ple,  the  tails  of  several  protein  absorption  bands  that 
were  not  included  in  the  modeling  might  contribute  to 
absorption  in  this  wavelength  region.  A  more  exten¬ 
sive  spectral  library  is  therefore  desirable  for  future 
studies.  Alternatively,  one  could  model  the  back¬ 
ground  empirically,  either  with  predetermined  math¬ 
ematical  functions  of  wavelen^h  (see,  e.g.,  Hull  et 
or  with  a  principal-component  analysis  of  sev¬ 
eral  background  residuals  once  many  samples  have 
been  studied.^® 

For  the  sake  of  comparison,  we  used  the  same 
least-squares  fitting  algorithm  to  calculate  the  spec¬ 
tral  fits  and  concentrations  twice,  once  with  the  full 
spectrum  and  once  with  only  the  FD  wavelengths, 
even  though  the  second  procedure  does  not  fully  ex¬ 
ploit  the  potential  of  the  SSFD  method.  Interesting 
differences  are  found  between  the  FD  and  SSFD  cal¬ 
culations.  Figure  7  and  Table  1  show  that  FD  sig¬ 
nificantly  underestimates  the  fat  concentration  in  the 
first  breast  sample:  0.75  g/cm^  with  the  SSFD  as 
opposed  to  0.44  g/cm^  with  the  FD.  The  reason  for 
this  discrepancy  lies  in  the  fact  that  no  laser  wave¬ 
length  is  close  enough  to  the  fat  peak  (at  present, 
commercial  laser  diodes  at  this  particular  wave¬ 
length  aue  uncommon).  Thus,  with  no  weighting  in 
this  region,  the  FD  fit  tolerates  large  fitting  errors 


near  the  absorption  maximum  of  fat  (as  illustrated  in 
Fig.  7),  leading  to  large  errors  in  the  fat  concentra¬ 
tion.  Such  a  discrepancy  is  naturally  enhanced  for 
this  sample,  for  which  the  fat  concentration  is  high¬ 
est.  This  example  illustrates  the  shortcomings  of 
using  a  limited  number  of  sources  (i.e.,  wavelengths), 
especially  when  a  significant  chromophore  lies  in  an 
undersampled  spectral  region.  A  similar  error  is 
seen  in  the  fitting  of  the  areolar  breast  location  (Fig. 
8);  this  time  it  is  the  water  peak  that  is  poorly  ad¬ 
dressed  by  the  FD  analysis. 

The  total  hemoglobin  determination  is  more  con¬ 
sistent  between  the  FD  and  SSFD  calculations,  owing 
to  the  smoother  spectral  features  of  these  chro¬ 
mophores  and  better  diode  coverage  of  the  650- 
850-nm  regime.  However,  the  oxygenation 
determination  by  FD  is  systematically  higher.  As 
discussed  above,  consistent  small  differences  occur 
between  the  fit  and  the  experiments  in  the  650- 
800-nm  region  because  of  some  incompleteness  in  our 
library  of  fitting  line  shapes.  Fitting  the  curve  with 
a  limited  number  of  wavelengths  naturally  exacer¬ 
bates  the  incompleteness  and  produces  a  bias  in  the 
results. 

As  mentioned  above,  Cubeddu  et  al.,^  working  in 
the  time  domain,  obtained  similar  optical  proper¬ 
ties  and  physiological  parameters  for  premeno¬ 
pausal  breast.  Conceptually,  the  two  experiments 
exploit  the  same  optical  phenomena  and  are  essen¬ 
tially  complementary  techniques  in  different  do¬ 
mains.  A  technical  advantage  of  the  SSFD 
measurement  is  the  use  of  a  white-light  source 
rather  than  a  tunable  laser.  Thus  the  SS  measure¬ 
ment  covers  all  wavelengths  simultaneously, 
whereas  the  laser(s)  must  be  tuned  separately  to 
each  wavelength.®  Additionally,  the  heart  of  the 
SSFD  system  (a  few  laser  diodes,  a  frequency  gen¬ 
erator,  an  avalanche  photodiode,  a  network  ana¬ 
lyzer  or  lock-in  amplifier,  and  a  SS  reflectance 
system)  is  inexpensive  and  easy  to  maintain  com¬ 
pared  with  a  tunable  laser  system  and  a  single- 
photon-counting  detection  apparatus. 

As  was  already  stated,  the  method  presented  here 
has  the  advantage  of  being  compatible  with  a  sin¬ 
gle,  large  source- detector  separation  (»10  trans¬ 
port  mean  free  paths).  In  contrast  to  spatially 
resolved  methods,  it  is  well  suited  for  interrogating 
deep  structures  in  relatively  heterogeneous  sam¬ 
ples.  As  numerous  authors  have  shown,!®  ®!  ®®  the 
layered  structure  of  tissue  affects  reflectance  differ¬ 
ently  at  different  source- detector  separations,  rais¬ 
ing  doubts  about  the  applicability  of  spatially 
resolved  techniques  that  assume  sample  homoge¬ 
neity  over  a  large  range.  Whereas  variations  in 
the  FD  modulation  frequency  w  do  change  the  op¬ 
tically  explored  tissue  volume,  these  effects  are 
modest  in  comparison  with  changes  in  p,  particu¬ 
larly  at  high  absorption.®!  Consequently,  the  as¬ 
sumption  of  homogeneity  is  less  extreme  for  the 
single-distance  measurements  proposed  here:  Its 
essence  is  that  all  measurements  at  one  wavelength 
explore  a  more  consistent  volume;  this  does  not 
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mean  that  different  wavelengths  will  explore  the 
same  volume.  Indeed,  less-absorbed  wavelengths 
will  explore  larger  regions  than  wavelengths  that 
attenuate  more  rapidly. 

7.  Conclusion 

A  combination  of  SS  and  FD  reflectance  measure¬ 
ments  has  been  described  for  absorption  spectroscopy 
of  turbid  media;  beneficial  aspects  of  both  techniques 
have  been  described.  As  with  SS,  the  wavelength 
coverage  is  continuous,  detecting  absorption  features 
that  may  not  be  discemable  in  the  discrete  wave¬ 
lengths  chosen  for  FD.  The  prediction  of  constituent 
concentrations,  for  instance,  in  breast  tissue,  is  sub¬ 
stantially  improved  when  full-spectrum  absorption 
data  rather  than  a  handful  of  wavelengths  are  used. 
As  with  FD,  however,  only  a  single  source- detector 
separation  is  required,  making  the  technique  more 
amenable  to  reporting  volume-averaged  values  for 
heterogeneous  samples.  In  addition,  the  source- 
detector  separation  can  be  large,  allowing  for 
centimeter-scale  mean  probing  depths  that  cannot  be 
achieved  with  spatially  resolved  SS  techniques. 
This  advantage  for  deeply  probing  studies  is  signifi¬ 
cant  for  many  clinical  purposes.  An  application  to 
breast  analysis  has  been  demonstrated,  with  quanti¬ 
tative  in  vivo  spectra  of  human  breast  obtained  rap¬ 
idly  (<  1  min) .  The  total  hemoglobin  content,  oxygen 
saturation,  and  water  and  fat  content  of  the  breast 
samples  have  been  calculated  from  the  spectra,  and 
failures  of  FD-only  fitting  have  been  highlighted. 
The  technique  is  relatively  inexpensive  and  could 
prove  valuable  for  improving  accuracy  in  the  devel¬ 
opment  of  quantitative  photon  migration  for  clinical 
instruments. 
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Il.d  Calculation  Tissue  Physiology:  Absorption  for  use  in  tissue  optics  studies,  although  it  did  not  change  the  essence  of  our  results  for  the 


IV.  DISCUSSION  IV.b  Normal  optical  properties  needed 
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